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2 DIFFERENTIALEQUATIONS 

Differential Equations 

A differential equation is an equation involving an unknown function and 
its derivatives. 

Example 1.1: The following are differential equations involving the un- 
known function y. 

Q- = 5x + ?> (1.1) 

dx 

e y^y. +2 (^] 2 =i (1.2) 

dx 2 \dx) 

4ll + ( s i nA ;)ll + 5^ = (1.3) 

dx dx 



y dx\ 



d y 



+ ^tx) +y \tx)= 5x (L4) 



dt 2 dx 



2 -4-f = (1.5) 



A differential equation is an ordinary differential equation if the unknown 
function depends on only one independent variable. If the unknown func- 
tion depends on two or more independent variables, the differential equa- 
tion is apartial differential equation. In this book we will be concerned 
solely with ordinary differential equations. 

Example 1.2: Equations 1.1 through 1.4 are examples of ordinary differ- 
ential equations, since the unknown function y depends solely on the vari- 
able x. Equation 1.5 is a partial differential equation, since y depends on 
both the independent variables t and x. 
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Note! 



The orderoí a differential equation is the order of 
the highest derivative appearing in the equation. 



Example 1.3: Equation 1.1 is a first-order differential equation; 1.2, 1.4, 
and 1.5 are second-order differential equations. (Note in 1.4 that the or- 
der of the highest derivative appearing in the equation is two.) Equation 
1.3 is a third-order differential equation. 

Notation 

The expressions y\y ,, ,y /// ,y (4 \...,y ( " ) are often used to represent, re- 
spectively, the first, second, third, fourth, . . ., nth derivatives of y with re- 
spect to the independent variable under consideration. Thus, y" repre- 

sents d 2 y I dx 2 if the independent variable is x, but represents d 2 y I dp 2 

if the independent variable is/?. Observe that parenthesis are used in y (n) 
to distinguish it from the nth power, y n . If the independent variable is 
time, usually denoted by t, primes are often replaced by dots. Thus, 

y, y, and y represent, dy/dt, d 2 y/dt 2 , and d 3 y/dt 3 , respectively. 

Solutions 

A solution of a differential equation in the unknown function y and the 
independent variable x on the interval 3 is a function y(x) that satisfies 
the differential equation identically for all x in 3 . 

Example 1.4: Is y(x) = c i sin2x + c 2 cos2x, where c, and c 2 are arbi- 
trary constants, a solution of y" + 4y = 0? 

Differentiating y, we find y' = 2c,cos2x - 2c 2 sin2x and y" = 
-4c , sin2x - 4c 2 eos 2a: . Henee, 
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y" + 4y = (-4c, sin2x-4c 2 cos2x) + 4(c, sin2.x: + c 2 cos2x) 
= (-4c, + 4c, ) sin 2x + (-4c 2 + 4c 2 ) eos 2x 
= 

Thus, y = c, sin 2x + c-, eos 2x satisfies the differential equation for all 
valúes of x and is a solution on the interval (-°o,°°). 

Example 1.5: Determine whether y = x 2 -1 is a solution of (y') 4 +y 2 
= -\. 

Note that the left side of the differential equation must be nonnegative for 
every real function y(x) and any x, since it is the sum of terms raised to 
the second and fourth powers, while the right side of the equation is neg- 
ative. Since no function y(x) will satisfy this equation, the given differ- 
ential equation has no solutions. 

We see that some differential equations have infinitely many solu- 
tions (Example 14), whereas other differential equations have no solu- 
tions (Example 1.5). It is also possible that a differential equation has ex- 
actly one solution. Consider (v') 4 +y 2 =0, which for reasons identical 
to those given in Example 1.5 has only one solution y = 0. 



You Need to Know 



/ 



A particular solution oí a differential equation is any 
one solution. The general solution of a differential 
equation is the set of all solutions. 



Example 1.6: The general solution to the differential equation in Ex- 
ample 1.4 can be shown to be (see Chapters Four and Five) v = 
c, sin2x + c 2 cos2x. That is, every particular solution of the differential 
equation has this general form. A few particular solutions are: (a) y = 
5sin2x-3cos2x (choose c, = 5 and c 2 =-3), (b) ;y = sin2x (choose 
c, = 1 and c 2 =0), and (c) y = (choose c ] =c 2 = 0). 
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The general solution of a differential equation cannot always be ex- 
pressed by a single formula. As an example consider the differential equa- 
tion y + y 2 = , which has two particular solutions y = 1 / x and y = . 

Initial-Value and Boundary-Value Problems 

A differential equation along with subsidiary 
conditions on the unknown function and its de- 
rivatives, all given at the same valué of the in- 
dependent variable, constitutes an initial-value 
problem. The subsidiary conditions are initial 
conditions. If the subsidiary conditions are giv- 
en at more than one valué of the independent 
variable, the problem is a boundary-value prob- 
lem and the conditions are boundary conditions. 

Example 1.7: The problem y" + 2v' = e x ;y{n) = \,y'{n) = 2 is an initial 
valué problem, because the two subsidiary conditions are both given at 
x = 7ü . The problem y" + 2 y' = e x ; y(0) = 1, y(\) = 1 is a boundary-value 
problem, because the two subsidiary conditions are given at x = and 
x=l. 

A solution to an initial-value or boundary-value problem is a func- 
tion y(x) that both solves the differential equation and satisfies all given 
subsidiary conditions. 

Standard and Differential Forms 

Standard form for a first-order differential equation in the unknown func- 
tion y(x) is 

/=/(*» y) (i.6) 

where the derivative y' appears only on the left side of 1.6. Many, but 
not all, first-order differential equations can be written in standard form 
by algebraically solving for y' and then setting f(x,y) equal to the right 
side of the resulting equation. 
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The right side of 1 .6 can always be written as a quotient of two oth- 
er functions M{x,y) and -N(x,y). Then 1.6 becomes dy I dx = M(x,y)/ 
—N(x,y), which is equivalent to the differential form 

M(x,y)dx + N(x,y)dy = (1.7) 

Linear Equations 

Consider a differential equation in standard form 1 .6. lff(x,y) can be writ- 
ten as f(x,y) = -p(x)y + q(x) (that is, as a function of x times y, plus an- 
other function of x), the differential equation is linear. First-order linear 
differential equations can always be expressed as 

y + p(x)y = q(x) (1.8) 

Linear equations are solved in Chapter Two. 

Bernoulli Equations 

A Bernoulli differential equation is an equation of the form 

y' + p(x)y = q(x)y n (1.9) 

where n denotes a real number. When n = 1 or n = 0, a Bernoulli equation 
reduces to a linear equation. Bernoulli equations are solved in Chapter 
Two. 

Homogeneous Equations 

A differential equation in standard form (1.6) is homogeneous if 

f(tx,ty) = f(x,y) (1.10) 

for every real number t. Homogeneous equations are solved in Chapter 
Two. 
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Note! 



In the general framework of differential equations, 
the word "homogeneous" has an entirely different 
meaning (see Chapter Four). Only in the context 
of first-order differential equations does "homoge- 
neous" have the meaning defined above. 



Separable Equations 

Consider a differential equation in differential form (1.7). If M(x,y) = A(x) 
(a function only of x) and N(x,y) = B(y) (a function only of y), the differ- 
ential equation is separable, or has its variables separated. Separable 
equations are solved in Chapter Two. 

Exact Equations 

A differential equation in differential form (1.7) is exact if 
dM(x,y) dN(x,y) 



dy dx 



(1.11) 



Exact equations are solved in Chapter Two (where a more precise defini- 
tion of exactness is given). 



Chapter 2 

Solutions of 
First-Order 

DlFFERENTIAL 

Equations 



In This Chapter: 

• Separable Equations 

• Homogeneous Equations 

• Exact Equations 

• Linear Equations 

• Bernoulli Equations 

• Solved Problems 

Separable Equations 
General Solution 

The solution to the first-order separable differential equation (see Chap- 
ter One). 

A(x)dx + B(y)dy = (2.1) 

8 
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is 

¡A(x)dx+JB(y)dy = c (2.2) 

where c represents an arbitrary constant. 

(See Problem 2.1) 

The integráis obtained in Equation 2.2 may be, for all practical pur- 
poses, impossible to evalúate. In such case, numerical techniques (see 
Chapter 14) are used to obtain an approximate solution. Even if the indi- 
cated integrations in 2.2 can be performed, it may not be algebraically 
possible to solve for y explicitly in terms of x. In that case, the solution is 
left in implicit form. 

Solutions to the Initial-Value Problem 

The solution to the initial- valué problem 

A(x)dx + B(y)dy = 0; y(x ) = y (2.3) 

can be obtained, as usual, by first using Equation 2.2 to solve the differ- 
ential equation and then applying the initial condition directly to evalú- 
ate c. 

Alternatively, the solution to Equation 2.3 can be obtained from 

x y 

\ A(s)ds + j B(t)dt = (2.4) 

*o yo 

where s and t are variables of integration. 

Homogeneous Equations 

The homogeneous differential equation 

^ = f(x,y) (2.5) 

dx 

having the property f(tx, ty) = f(x, y) (see Chapter One) can be trans- 
formed into a separable equation by making the substitution 
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y = xv (2.6) 

along with its corresponding derivative 

d y , clv (r\ n\ 

~¿- = v + x — (2.7) 

dx dx 

The resulting equation in the variables v and x is solved as a separable 
differential equation; the required solution to Equation 2.5 is obtained by 
back substitution. 

Alternatively, the solution to 2.5 can be obtained by rewriting the dif- 
ferential equation as 

dx 1 

(2.8) 



and then substituting 



dy f(x,y) 



x = yu (2.9) 



and the corresponding derivative 

dx du 

— = u + y — (2.10) 

dy dy 

into Equation 2.8. After simplifying, the resulting differential equation 
will be one with variables (this time, u and y) separable. 

Ordinarily, it is immaterial which method of solution is used. Occa- 
sionally, however, one of the substitutions 2.6 or 2.9 is definitely superi- 
or to the other one. In such cases, the better substitution is usually appar- 
ent from the form of the differential equation itself. 

(See Problem 2.2) 

Exact Equations 

Defining Properties 

A differential equation 

M(x, y)dx + N(x, y)dy = (2.11) 
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is exact if there exists a function g(x, y) such that 

dg(x, y) = M(x, y)dx + N(x, y)dy (2.12) 



^ 



Note! 



Test for exactness: If M(x,y) and N(x,y) are con- 
tinuous functions and have continuous first partial 
derivatives on some rectangle of the xy-plane, then 
Equation 2.11 is exact if and only if 
dM(x,y) dN(x,y) 



dy dx 



(2.13) 



Method of Solution 

To solve Equation 2.11, assuming that it is exact, first solve the equations 

dg(x,y) 



dx 

dg(*,y) 



= M(x,y) (2.14) 

= N(x,y) (2.15) 



dy 

for g(x, y). The solution to 2. 1 1 is then given implicitly by 

g(x,y) = c (2.16) 

where c represents an arbitrary constant. 

Equation 2. 16 is immediate from Equations 2. 1 1 and 2. 12. If 2. 12 is 
substituted into 2. 1 1 , we obtain dg(x, y(x)) = 0. Integrating this equation 
(note that we can write as dx), we have \ dg(x,y(x)) = \0dx, which, 
in turn, implies 2.16. 

Integrating Factors 

In general, Equation 2. 1 1 is not exact. Occasionally, it is possible to trans- 
form 2. 1 1 into an exact differential equation by a judicious multiplica- 
tion. A function I(x, y) is an integrating factor for 2. 1 1 if the equation 
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I(x, y) [M(x, y)dx + N(x, y)dy] = (2.17) 

is exact. A solution to 2. 1 1 is obtained by solving the exact differential 
equation defined by 2.17. Some of the more common integrating factors 
are displayed in Table 2. 1 and the conditions that follow: 



ifl 

N 



r 'dM dN 



= g(x), a function of x alone, then 
dy dx 



I(x,y) = J g{x)(lx (2.18) 



Ifi- 

M 



r óM dN 



= h(y), a function ofy alone, then 
dy dx 



I{x,y) = e^ h(y)dy (2.19) 

If M = yf(xy) and N = xg(xy), then 

I(x,y) = ! (2.20) 

xM-yN 

In general, integrating factors are difficult to unco ver. If a differential 
equation does not have one of the forms given above, then a search for 
an integrating factor likely will not be successful, and other methods of 
solution are recommended. 

(See Problems 2.3-2.6) 

Linear Equations 

Method of Solution 

A first-order linear differential equation has the form (see Chapter One) 

y' + p(x)y = q(x) (2.21) 

An integrating factor for Equation 2.21 is 

I(x) = e J (2.22) 
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Table 2.1 



Group of terms 



Integrating factor l(x, y) 



Exact differential dy{x,y) 



y dx - x dy 



x dy —ydx_ I y 
x 2 \x 



y dx — x dy 



y dx — x dy _ fx 



r 



y dx — x dy 



xy 



xdv-vdx ,/, v 
— : = d ln - 

XV V X 



v dx - xdv 



* 2 + r 



xdy-ydx í y 

— ^ — -^— = d[ arctan - 

x~ + y \ xí 



y dx + x dy 



1 
xv 



V dx + x dy 
xy 



= d(lnxy) 



y dx 4- x dy 



(xy)'" 



n>l 



y dx + x dy _ I" — 1 

~~(xjy~ = C [{^-\)(¿y~)"~. 



v dy + x dx 



vdy +x dx ri 

x +y 12 



y dy + x dx 



(jr+y 2 )"' 



n>\ 



y dy + x dx í" — 1 

7^+77" ~ d l2ln-\)(x T +y r f 



ay dx + bx dy 
(a, b constants) 



x" y 



x" y \aydx+bxdy) = d(x"y") 



which depends only onx and is independent ofy. When both sides of 2.21 
are multiplied by I(x), the resulting equation 



I(x)y' + p(x)I(x)y = I(x)q(x) 



(2.23) 



is exact. This equation can be solved by the method described previous- 
ly. A simpler procedure is to rewrite 2.23 as 



d(yl) 
dx 



= Iq(x) 
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intégrate both sides of this last equation with respect to x, and then solve 
the resulting equation for y. The general solution for Equation 2.21 is 

í I(x)q(x)dx + c 

y = - 

I(x) 

where c is the constant of integration. 
(See Problem 2.7) 

Bernoulli Equations 

A Bernoulli differential equation has the form 

y' + p(x)y = q(x)y n (2.24) 

where n is a real number. The substitution 

z = y l ~ n (2.25) 

transforms 2.24 into a linear differential equation in the unknown func- 
tion z(x). 

(See Problem 2.8) 

Solved Problems 

Solved Problem 2.1 Solve ^- = X + . 

dx y 

This equation may be rewritten in the differential form 
(x 2 +2)dx-ydy = 
which is separable with A (x) = x 2 + 2 and B(y) = -y. Its solution is 
\(x 2 +2)dx- \ydy = c 



or 



-x 3 +2x--y 2 =c 

3 2 J 
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Solving for y, we obtain the solution in implicit form as 

2 
y 2 =— x 3 +4x + k 
y 3 

with k = -2c. Solving for y explicitly, we obtain the two solutions 



y= l—x 3 +4x + k and y = - I— x 3 +4x + k 

Solved Problem 2.2 Solve / = 2±£. 

x 

This differential equation is not separable. Instead it has the form 
y' = f(x,y), with 

fí \ y +x 

f(x,y) = 

x 

where 

ty + tx t(y + x) y + x 
f(tx,ty) = 2. = M i = l = f( x ,y) 

tx tx X 

so it is homogeneous. Substituting Equations 2.6 and 2.7 into the equa- 
tion, we obtain 

dv xv + x 

V + x — = 

dx x 

which can be algebraically simplified to 

dv 1 

x — = 1 or — dx — dv = 
dx x 

This last equation is separable; its solution is 

| — dx — \ dv = c 
which, when evaluated, yields v = ln I x I — c, or 

v = \n\kx\ (2.26) 

where we have set c = -ln I A: I and have noted that ln \x l + ln I £ 1 = ln I xk I. 
Finally, substituting v = y/x back into 2.26, we obtain the solution to the 
given differential equation as y = x ln I kx I . 
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Solved Problem 2.3 Solve Ixydx + (1 + x 2 )dy = 0. 

This equation has the form of Equation 2. 1 1 with M(x, y) = 2xy and 
N(x, y) = 1 + x 2 . Since dM I dy = 6N I dx = 2x, the differential equation 
is exact. Because this equation is exact, we now determine a function 
g(x, y) that satisfies Equations 2.14 and 2.15. Substituting M(x, y) = 2xy 
into 2. 14, we obtain óg / dx = 2xy. Integrating both sides of this equation 
with respect to x, we find 



[JLdx=¡ Ixydx 

J rtX J 



or 



g(x,y) = x 2 y + h(y) (2.27) 

Note that when integrating with respect to x, the constant (with respect to 
x) of integration can depend on y. 

We now determine h(y). Differentiating 2.27 with respect to y, we 
obtain óg / dy = x 2 + h'(y) Substituting this equation along with N(x, y) = 
1 + x 2 into 2. 15, we have 

x 2 +h\y) = \ + x 2 or tí(y) = \ 

Integrating this last equation with respect to y, we obtain h(y) =y + c 1 
(c { = constant). Substituting this expression into 2.27 yields 

g(x,y) = x 2 y + y + c { 

The solution to the differential equation, which is given implicitly by 
2.16 as g(x, y) = c, is 

x 2 y + y = c 2 (c 2 = c-c x ) 

Solving for y explicitly, we obtain the solution as y = c 2 /(x 2 +1). 

Solved Problem 2.4 Determine whether the differential equation ydx - 
xdy = is exact. 
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This equation has the form of Equation 2. 1 1 with M (x, y) = y and N (x, y) 
= -x. Here 

dM . . dN 

= 1 and — = -1 

dy dx 

which are not equal, so the differential equation is not exact. 

Solved Problem 2.5 Determine whether -\/x 2 is an integrating factor 
for the differential equation ydx - xdy = 0. 

It was shown in Problem 2.4 that the differential equation is not exact. 
Multiplying it by -l/x 2 , we obtain 

^ T (ydx-xdy) = or ^-dx + -dy = (2.28) 

X X X 

Equation 2.28 has the form of Equation 2.11 with M (x, y) = —y/x 2 and 
N (x, y) = l/x. Now 

dM _ 3f-/L-l_ 3 ÍO_ dN 

dy óy\x 2 ) x 2 dx\x) dx 

so 2.28 is exact, which implies that -\/x 2 is an integrating factor for the 
original differential equation. 

Solved Problem 2.6 Solve ydx - xdy = 0. 

Using the results of Problem 2.5, we can rewrite the given differential 
equation as 

xdy - ydx _ 

~V 2 

which is exact. Equation 2.28 can be solved using the steps described in 
Equations 2.14 through 2.16. 

Alternatively, we note from Table 2.1 that 2.28 can be rewritten as 
d (y/x) = 0. Henee, by direct integration, we have y/x = c, or y = xc, as 
the solution. 

Solved Problem 2.7 Solve y' + (A/x)y = x 4 . 
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The differential equation has the form of Equation 2.21, with /?(.*) = Al x 
and q(x) = x 4 , and is linear. Here 

jp(x)dx = í -dx = 41n I x I = ln x 4 

so 2.22 becomes 

I(x) = J p(x)c ' x =e lnx4 =x 4 (2.29) 

Multiplying the differential equation by the integrating factor de- 
fined by 2.29, we obtain 

x 4 y' + 4x 3 y = x 8 or — (yx 4 ) = x* 
dx 

Integrating both sides of this last equation with respect to x, we obtain 



4 1 o C 1 5 

= — x +c or v = — r + — x 
9 y x 4 9 



Solved Problem 2.8 Solve y' + xy = xy 2 . 

This equation is not linear. It is, however, a Bernoulli differential equa- 
tion having the form of Equation 2.24 with p(x) = q{x) = x, and n = 2. We 
make the substitution suggested by 2.25, namely z = >' 1-2 = > ,_1 , from 
which follow 

1 A > Z ' 

y = — and y = — - 

z z~ 

Substituting these equations into the differential equation, we obtain 

z' x x 
— - + — = ^r or z -xz = -x 

z z z 

This last equation is linear for the unknown function z(x). It has the form 
of Equation 2.21 with y replaced by z and p(x) = q(x) = -x. The integrat- 
ing factor is 

I {x) = J ( - x)dx =e- x2 ' 2 
Multiplying the differential equation by I{x), we obtain 
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e -* 2 <2Ck_ xe - X > / 2 z = _ xe - X >,2 

dx 



or 



dx\ I 

Upon integrating both sides of this last equation, we have 

ze =e +c 



whereupon 



Z(x) = ce x /2 +l 



The solution of the original differential equation is then 

_J__ 1 



Z ce x/2 +l 
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Growth and Decay Problems 

Let N(t) denote the amount of substance (or popula- 
tion) that is either growing or decaying. If we assume 
that dNIdt, the time rate of change of this amount of 
substance, is proportional to the amount of substance 
present, then dNIdt = kN, or 
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— -kN = (3.1) 

dt 

where k is the constant of proportionality. 

We are assuming that N(t) is a differentiable, henee continuous, func- 
tion of time. For population problems, where N(t) is actually discrete and 
integer-valued, this assumption is incorrect. Nonetheless, 3.1 still pro- 
vides a good approximation to the physical laws governing such a sys- 
tem. 

Temperature Problems 

Newton's law of cooling, which is equally applicable to heating, states 
that the time rate ofehange ofthe temperature of a body is proportional 
to the temperature difference between the body and its surrounding médi- 
um. Let T denote the temperature of the body and let T m denote the tem- 
perature of the surrounding médium. Then the time rate of change of the 
temperature ofthe body is dT/dt, and Newton's law of cooling can be for- 
mulated as dT/dt = —k(T -T ), or as 

^L + kT = kT m (3.2) 

dt 

where k is a positive constant of proportionality. Once k is chosen posi- 
tive, the minus sign is required in Newton's law to make dT/dt negative 
in a cooling process, when T is greater than T m , and positive in a heating 
process, when T is less than T . 

Falling Body Problems 

Consider a vertically falling body of mass m that is being influenced only 
by gravity g and an air resistance that is proportional to the velocity of 
the body. Assume that both gravity and mass remain constant and, for 
convenience, choose the downward direction as the positive direction. 
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You Need to Know 



/ 



Newton 's second law of motion: The net forcé 
acting on a body is equal to the time rate ofchange 
of the momentum of the body; or, for constant 
mass, 

r- dV 

F=m - Q J ( 33 ) 

where F is the net forcé on the body and v is the 
velocity of the body, both at time t. 



For the problem at hand, there are two forces acting on the body: the forcé 
due to gravity given by the weight w of the body, which equals mg, and 
the forcé due to air resistance given by —kv, where k > is a constant of 
proportionality. The minus sign is required because this forcé opposes the 
velocity; that is, it acts in the upward, or negative, direction (see Figure 
3-1). The net forcé F on the body is, therefore, F = mg — kv. Substituting 
this result into 3.3, we obtain 

mg-kv = m — 
5 dt 

or 

dv k ,„ „ 

— + -v = g (3.4) 

dt m 

as the equation of motion for the body. 

If air resistance is negligible or nonexistent, then k = and 3.4 sim- 
plifies to 

When k > 0, the limiting velocity v / is defined by 

v,=^ (3.6) 
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Falling body 




, ' Positive x-direct¡on 
Figure 3-1 



Caution: Equations 3.4, 3.5, and 3.6 are valid only if the given con- 
ditions are satisfied. These equations are not valid if, for example, air re- 
sistance is not proportional to velocity but to the velocity squared, or if 
the upward direction is taken to be the positive direction. 

Dilution Problems 



Consider a tank which initially holds V Q gal of brine that contains a Ib of 
salt. Another solution, containing b Ib of salt per gallón, is poured into the 
tank at the rate of e gal/min while simultaneously, the well-stirred solu- 
tion leaves the tank at the rate of/ gal/min (Figure 3-2). The problem is 
to find the amount of salt in the tank at any time t. 

Let Q denote the amount (in pounds) of salt in the 
tank at any time. The time rate of change of Q, dQ/dt, 
equals the rate at which salt enters the tank minus the rate 
at which salt leaves the tank. Salt enters the tank at the 
rate of be Ib/min. To determine the rate at which salt 
leaves the tank, we first calcúlate the volume of brine 
in the tank at any time t, which is the initial volume V Q 
plus the volume of brine added et minus the volume of brine removed ft. 
Thus, the volume of brine at any time is 
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V +et-ft 



(3.7) 



The concentration of salt in the tank at any time is Q I ( V +et — ft), from 
which it follows that salt leaves the tank at the rate of 



Thus, 



or 



( 



f 



c}Q 
dt 



Q 



Vo+et-ft 



= be-f 



lb/min 



Q 



V +et-ft 



f 



dt V +(e-f)t 



Q = be 



(3.8) 



e gal/min 




/gal/min 



Figure 3-2 
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Figure 3-3 



Electrical Circuits 



The basic equation governing the amount of current / (in amperes) in a 
simple RL circuit (see Figure 3-3) consisting of a resistance R (in ohms), 
an inductor L (in henries), and an electromotive forcé (abbreviated emf ) 
E (in volts) is 



di R T E 
— + — / = — 
dt L L 



(3.9) 



For an RC circuit consisting of a resistance, a capacitance C (in farads), 
an emf, and no inductance (Figure 3-4), the equation governing the 
amount of electrical charge q (in coulombs) on the capacitor is 



da 1 E 

-i + q = — 

dt RC R 



(3.10) 



R 



Figure 3-4 
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The relationship between q and / is 

I= d -Í (3.11) 

dt 

For more complex circuits see Chapter Seven. 

Orthogonal Trajectories 

Consider a one-parameter family of curves in the xy-plane defined by 

F(x,y,c) = (3.12) 

where c denotes the parameter. The problem is to find another one-pa- 
rameter family of curves, called the orthogonal trajectories of the fami- 
ly of curves in 3.12 and given analytically by 

G(x,y,k) = (3.13) 

such that every curve in this new family 3.13 intersects at right angles 
every curve in the original family 3.12. 

We first implicitly differentiate 3. 12 with respect to x, then eliminate 
c between this derived equation and 3.12. This gives an equation con- 
necting x, y, and y', which we solve for y' to obtain a differential equa- 
tion of the form 

T^ = /(*.?) (3.14) 

dx 



The orthogonal trajectories of 3.12 are the solutions of 

dy 1 



(3.15) 



dx f(x,y) 

For many families of curves, one cannot explicitly solve for dyldx 
and obtain a differential equation of the form 3.14. We do not consider 
such curves in this book. 

Solved Problems 

Solved Problem 3.1 A bacteria culture is known to grow at a rate pro- 
portional to the amount present. After one hour, 1 000 strands of the bac- 
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teria are observed in the culture; and after four hours, 3000 strands. Find 
(a) an expression for the approximate number of strands of the bacteria 
present in the culture at any time t and (b) the approximate number of 
strands of the bacteria originally in the culture. 

(a) Let N(t) denote the number of bacteria strands in the culture at 
time t. From Equation 3. 1 , dNIdt - kN = 0, which is both linear and sep- 
arable. Its solution is 



At t = 1, N = 1000; henee, 



At t = 4,N = 3000; henee. 



N(t) = ce kt (3.16) 



1000 = ce k (3.17) 



3000 = ce 4 * (3.18) 



Solving 3.17 and 3.18 for k and c, we find 

k = -\n3 « 0.3662 and c = 1000^ = 693 
3 

Substituting these valúes of A: and c into 3.16, we obtain 

N(t) = 693e 03662T (3.19) 

(b) We require Natt = 0. Substituting t = into 3. 19, we obtain N(0) 
= 693é>< a3662 x°> = 693. 

Solved Problem 3.2 A tank initially holds 100 gal of a brine solution 
containing 20 Ib of salt. At t = 0, fresh water is poured into the tank at the 
rate of 5 gal/min, while the well-stirred mixture leaves the tank at the 
same rate. Find the amount of salt in the tank at any time t. 

Here, V Q = 100, a = 20, b = 0, and e =f= 5. Equation 3.8 becomes 
dt 20 
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The solution of this linear equation is 

Q = ce'" 20 (3.20) 

At t = 0, we are given that Q = a = 20. Substituting these valúes into 3.20, 
we find that c = 20, so that 3.20 can be rewritten as Q = 20e~ t/2 °. Note 
that as t — > «o, Q — > as it should, since only fresh water is being added. 
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• Linear Differential Equations 
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• The Wronskian 

• Nonhomogeneous Equations 

Linear Differential Equations 

An nth-order linear differential equation has the form 

b n (x)y (n) +b n _ i (x)y ( "- i) + - + b l (x)y' + b (x)y = g(x) (4.1) 

where g(x) and the coefficients b .(x) (j = 0,1 ,2,..., ri) de- 
pend solely on the variable x. In other words, they do 
not depend on y or any derivative of y. 
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If g(x) = 0, then Equation 4. 1 is homogeneous; if not, 4. 1 is nonho- 
mogeneous. A linear differential equation has constant coefficients if all 
the coefficients b (x) in 4. 1 are constants; if one or more of these coeffi- 
cients is not constant, 4. 1 has variable coefficients. 

Theorem 4.1. Consider the initial- valué problem given by the linear dif- 
ferential equation 4. 1 and the n initial conditions 



y(x ) = c , y'(x ) = c l , 
f{x ) = c 2 ,...,y^\x Q ) = c n _ x 



\f g{x) and b(x) (j = 0,1,2,..., n) are continuous in some interval $ con- 
taining x Q and if b n {x) ■£ in $, then the initial- valué problem given by 
4.1 and 4.2 has a unique (only one) solution defined throughout $. 

When the conditions on b (x) in Theorem 4. 1 hold, we can divide 
Equation 4. 1 by b n {x) to get 

y (n) + a n _ x (x)y (n - l) +- + a 2 (x)f + a, (x)y' + a (x)y = <¡>(x) (4.3) 

where aXx) = b(x)/bjx) (j = 0,1,2,..., n - 1) and <¡)(x) = g(x)/b n (x). 
Let us define the differential operator h(y) by 

L( v) ^ y (n) + a n _ x (x) v ( "-' > + • • • + a 2 {x)f + a x {x)y' + a (x)y (4.4) 

where a.(x) (i = 0,1,2,..., n — 1) is continuous on some interval of interest. 
Then 4.3 can be rewritten as 

L();) = 0(x) (4.5) 

and, in particular, a linear homogeneous differential equation can be ex- 
pressed as 

L0;) = (4.6) 
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Linearly Independent Solutions 

A set of functions {y l (x),y 2 (x),...,y n (x)} is linearly clependent on a <x<b 
if there exist constants c l ,c 2 ,...,c n not all zero, such that 



c 1 y 1 (x) + c 2 y 2 (x) + --- + c n y n (x) = 



(4.7) 



on a <x < b. 



Example 4.1: The set {x,5x,l,sin x}is linearly dependent on [-1,1] since 
there exist constants c y = -5, c 2 = 1 , c 3 = 0, and c 4 = 0, not all zero, such 
that 4.7 is satisfied. In particular, 

-5 • x + 1 • 5x + • 1 + • sin x = 

Note that c , = c 2 = • • • c n = is a set of constants that always satis- 
fies 4.7. A set of functions is linearly dependent if there exists another set 
of constants, not all zero, that also satisfies 4.7. If the only solution to 4.7 
is c, = c 2 = ■■■ c n = 0, then the set of functions {yfa), y 2 (x),...,y n (x)} is 
linearly independent on a < x < b. 

Theorem 4.2. The nth-order linear homogeneous differential equation 
L(y) = always has n linearly independent solutions. ]fy 1 (x),y 2 (x),...,y n (x) 
represent these solutions, then the general solution of L(;y) = is 

y(x) = cj^x) + c 2 y 2 (x) +■ ■ ■+ cj n (x) (4.8) 

where c v c 2 ,...,c n denote arbitrary constants. 

The Wronskian 

The Wronskian of a set of functions {z^x), z 2 (x),..., Z n (x)} on the interval 
a<x < b, having the property that each function possesses n — 1 deriva- 
tives on this interval, is the determinant 



w{z x ,z 2 ,...,z n ) = 



zl 
z'l 



z 2 
zí 
zS 



(n-\) (»-l) 

C\ z 2 



Z n 
Z' n 

zl 

T ("-l) 
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Theorem 4.3. If the Wronskian of a set of n functions defined on the in- 
terval a < x < b is nonzero for at least one point in this interval, then the 
set of functions is linearly independent there. If the Wronskian is identi- 
cally zero on this interval and if each of the functions is a solution to the 
same linear differential equation, then the set of functions is linearly de- 
pendent. 

Caution: Theorem 4.3 is silent when the Wronskian is identically 
zero and the functions are not known to be solutions of the same linear 
differential equation. In this case, one must test directly whether Equa- 
tion 4.7 is satisfied. 

Nonhomogeneous Equations 

Let v denote any particular solution of Equation 4.5 (see Chapter One) 
and let y h (henceforth called the homo gene ous or complementan/ solu- 
tion) represent the general solution of the associated homogeneous equa- 
tion L(v) = 0. 

Theorem 4.4. The general solution to L(v) = (¡){x) is 



y=y h + y a 



(4.9) 



Don't Forget 

The general solution is the sum of 
the homogeneous and particular so- 
lutions. 
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The Characteristic Equation 

Second-Order Equations 

Corresponding to the differential equation 

y" + a iy ' + a o y = (5.1) 

in which a i and a Q are constants, is the algebraic equation 

X 2 +a l X + a o =0 (5.2) 

which is obtained by substituting y = e** (assuming x to be the indepen- 
dent variable) into Equation 5.1 and simplifying. Note that Equation 5.2 
can be easily obtained by replacing y", y', and y by X 2 , V, and, X o = 1, 
respectively. Equation 5.2 is called the characteristic equation of 5.1. 

Example 5.1 The characteristic equation of y" + 3y' - 4y = is X 2 + 31 
-4 = 0; the characteristic equation of y" - 2y' + y = is X 2 - 21 + 1 = 0. 

Characteristic equations for differential equations having dependent 
variables other than y are obtained analogously, by replacing the y'th de- 
rivan ve of the dependent variable by V (j = 0,1 ,2). 

nth-Order Equations 

Similarly, the characteristic equation of the differential equation 

y (n) +<W ( " _1) +- + a i y + a o y = (5.3) 

with constant coefficients a.(j = 0,1,..., n— 1) is 

¿" + a n _ i r- ] +--- + a x ¿ + a o =0 (5.4) 

The characteristic equation 5.4 is obtained from 5.3 by replacing y^ by 
V {j = 0,1,..., n). Characteristic equations for differential equations hav- 
ing dependent variables other than y are obtained analogously, by 
replacing they'th derivative of the dependent variable by AJ (j = 0,1,..., 
n). 
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Example 5.2. The characteristic equation of )> (4) - 3/" + 2y" — y = is 
X a — 3X 3 + 2X 2 -1=0. The characteristic equation of 

d 5 x d 3 x dx 

— s--3 — r + 5 lx = \) 

dr dt dt 



is 



¿ 5 -3A 3 +52-7 = 



# 



Caution! 



Characteristic equations are only defined for linear 
homogeneous differential equations with constant 
coefficients. 



General Solution for Second-Order Equations 

The characteristic equation (5.2) can be factored into 

(2-/i 1 )(/i-2 2 ) = (5.5) 



i£ 



Note! 



The roots of the characteristic polynomial deter- 
mine the solution of the differential equation. 



There are three cases to consider. 

Case 1. X t and X 2 both real and distinct. Two linearly indepen 
dent solí 
rem 4.2) 



dent solutions are e x,x and e x - x , and the general solution is (Theo- 
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y = c i e iX +c 2 e (5.6) 

In the special case X 2 = -X x , the solution 5.6 can be rewritten as y = 
k x cosh X x x + k 2 sinh X x x. 

Case 2. X í =a + ib,a complex number. Since a x and a Q in 5. 1 and 
5.2 are assumed real, the roots of 5.2 must appear in conjúgate pairs; 
thus, the other root is X 2 = a - ib. Two linearly independent solutions 
are e (a+l )x , and e (a ~' )x , and the general complex solution is 

y = d/ a+ib)x +d 2 e ( °- ib)x (5.7) 

which is algebraically equivalent to 

y = c 1 e ax cosbx + ^e^sinbx (5.8) 

Case 3. X 1 = X 2 . Two linearly independent solutions are e ,x and 
xe^ ,x , and the general solution is 



y = c x e Z]X +c 2 xe X]X (5.9) 

Warning: The above solutions are not valid if the differential equa- 
tion is not linear or does not have constant coefficients. Consider, for ex- 
ample, the equation y" - x 2 y = 0. The roots of the characteristic equa- 
tion are Aj = x and X 2 = -x, but the solution is not 

y = c / x)x + c 2 e ( ~ x)x = c x e x2 + c 2 e~ x2 

Linear equations with variable coefficients are considered in Chapter 
Twelve. 



General Solution for nth-Order Equations 

The general solution of 5.3 is obtained directly from the roots of 5.4. If 
the roots Aj, A 2 ,..., A^are all real and distinct, the solution is 

y = c x e x * x +c 2 e x ' x +--- + c n e x " x (5.10) 




CHAPTER 5: Linear Homogeneous Differential Equations 37 

If the roots are distinct, but some are com- 
plex, then the solution is again given by 5.10. 
As in the second-order equation, those terms 
involving complex exponentials can be com- 
bined to yield terms involving sines and co- 
sines. If X k is a root of multiplicity p [that is, if 
(X - X k ) p is a factor of the characteristic equa- 
tion, but (X - X k ) ,y+X is not] then there will be 
p linearly independent solutions associated 
with X k given by e kX ,xe kX ,x 2 e k * ,...,x p ~ x e kX . These solutions are 
combined in the usual way with the solutions associated with the other 
roots to obtain the complete solution. 

In theory it is always possible to factor the characteristic equation, 
but in practice this can be extremely difficult, especially for differential 
equations of high order. In such cases, one must often use numerical tech- 
niques to approximate the solutions. See Chapter Fourteen. 

Solved Problems 

Solved Problem 5.1 Solve y"-y'-2y = 0. 

The characteristic equation is X 2 - X - 2 = 0, which can be factored 
into (X -2)(X + 1) = 0. Since the roots X x =2 and X = -1 are real and dis- 
tinct, the solution is given by 5.6 as 

y = c x e lx + c 2 e~ x 

Solved Problem 5.2 Solve y"-8/ + 16y = 0. 

The characteristic equation is X 2 -82 + 16 = which can be fac- 
tored into (X - 4) 2 = 0. The roots X l = X 2 = 4 are real and equal, so the gen- 
eral solution is given by 5.9 as 

y = c ] e 4x +c 2 xe 4x 

Solved Problem 5.3 Solve /" - 6y" + 2/ + 36v = 0. 

The characteristic equation X 3 - 6X 2 + 2X + 36 = 0, has roots X x = 
—2, X 2 =4 + is¡2 , and X 3 =4 — /-y/2 . The solution is 
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y = c x e~ 2x + d 2 é 4+i42)x + d 3 e {4 -' ,42)x 
which can be rewritten, using Euler's relations 

e' bx = eos bx + i sin bx and e~ ibx = eos bx - i sin bx 
as 



y = c,e~ 2x + d 2 e 4x e H2x + d 3 e 4x e~ 

y = c x e~ 2x + d 2 e 4x (cos^x + /sin V2x) + d 3 e 4x (cos-j2x - isin42x) 

y = c x e~ 2x + (d 2 + d 3 )e 4x eos V2a: + i(d 2 - d 3 )e 4x sin 4lx 

y = c l e~ 2x +c 2 e 4x cos^¡2x + c 3 e 4x sinV2x 

Note that this form of the solution corresponding to the complex roots 
can be easily formulated using Equation 5.8. 
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The general solution to the linear differential equa- 

tion L(v) = (¡){x) is gi ven by Theorem 4.4 as y = y h + 

v where y denotes one solution to the differential 
p.p. 
equation and y h is the general solution to the asso- 

ciated homogeneous equation, h(y) = 0. Methods 
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for obtaining y h when the differential equation has constant coefficients 
are given in Chapter Five. In this chapter, we give methods for obtaining 
a particular solution y once y h is known. 

The Method of Undetermined Coefficients 

Simple Form of the Method 

The method of undetermined coefficients is applicable only if (p(x) and all 
of its derivatives can be written in terms of the same finite set of linearly 
independent functions, which we denote by {y^x), y 2 (x),..., y n (x)}. The 
method is initiated by assuming a particular solution of the form 

y p (x) = A 1 y 1 (x)+A 2 y 2 (x) + — + A n y n (x) 

where A , , A 2 ,...,A n denote arbitrary multiplicative constants. These arbi- 
trary constants are then evaluated by substituting the proposed solution 
into the given differential equation and equating the coefficients of like 
terms. 

Case 1. <¡>{x) =p n (x), an wth-degree polynomial in*. Assume a so- 
lution of the form 

y p =e ax {A n x n +A n _ x x"- x +- + A x x- (6.1) 

where A . (j = 0,1,2,..., n) is a constant to be determined. 

Case 2. (¡>{x) = he*** where k and a are known constants. Assume 

a solution of the form 

y p = Ae ax (6.2) 

where A is a constant to be determined. 

Case 3. Q(x) =k x sin fix + k 2 eos fix where k v k 2 , and fi are known 

constants. Assume a solution of the form 

y = A sin fix + B eos fix (6.3) 

where A and B are constants to be determined. 
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y = A sin px + B eos px in its entirety is assumed 
for 0(x) = /c, sin px+ k 2 eos px even when /t, or k 2 
is zero, because the derivatives of sines or cosines 
involve both sines and cosines. 



Generalizations 

If (¡){x) is the product of terms considered in Cases 1 through 3, take y to 
be the product of the corresponding assumed solutions and algebraically 
combine arbitrary constants where possible. In particular, if (j){x) = e™ pjx) 
is the product of a polynomial with an exponential, assume 



v =e ax (A n x n +A n _ l x"- l +- + A l x + A ) 



(6.4) 



where A . is as in Case 1 . If, instead, (¡){x) = e"* p n (x) sin fix is the product 
of a polynomial, exponential, and sine term, or if (¡)(x) = e ax p n (x) eos @x 
is the product of a polynomial, exponential, and cosine term, then assume 

y p = e ax sin/3x(A n x n +A n _ i x n - l +- + A i x + A ) 
+e ax cospx(B n x" + B n _ x x"~ l + ■ ■ ■ + £, x + B ) 



where A . and B(j = 0,1,2,..., n) are constants which still must be deter- 
mined. 

If 0(x) is the sum (or difference) of terms already considered, then 
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we take y to be the sum (or difference) of the corresponding assumed So- 
lutions and algebraically combine arbitrary constants where possible. 

Modifications 

If any term of the assumed solution, disregarding multiplicative con- 
stants, is also a term ofy h (the homogeneous solution), then the assumed 
solution must be modified by multiplying it by x m , where m is the small- 
est positive integer such that the product of x m with the assumed solution 
has no terms in common with y h . 

Limitations of the Method 

In general, if (¡)(x) is not one of the types of functions considered abo ve, 
or if the differential equation does not have constant coefficients, then the 
following method is preferred. 

Variation of Parameters 

Variation of parameters is another method for finding a particular solu- 
tion of the nth-order linear differential equation 

Uy) = *x) (6.6) 

once the solution of the associated homogeneous equation L(y) = is 
known. Recall from Theorem 4.2 that if y x (x), y 2 (x),..., y n (x) are n linear- 
ly independent solutions of L(y) = 0, then the general solution of L(y) = 
Oís 

>'ft= C iyiW + C 2) , 2W + - + C «)'«W ( 6 - 7 ) 



The Method 

A particular solution of L()>) = 0(x) has the form 

y p = v l y l +v 2 y 2 +--- + v n y n (6.8) 

where y i = yJ(x) (i = 1,2,..., n) is given in Equation 6.7 and v i (i = 1,2,..., 
n) is an unknown function of x which still must be determined. 
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To find v p first solve the following linear equations simultaneously 
for v-: 



V í>'l+ V 2) ? 2 +~' + v'„yn =0 

v\y\ + v 2 y 2 +- + v' n y' n = 

^r 2) +v^- 2 >+... + v^- 2 >=o 

v¡y¡ n - l) + v , 2 y ( 2- l) + -+v n y^=<i>(x) 



(6.9) 



Then intégrate each v\ to obtain v f , disregarding all constants of integra- 
tion. This is permissible because we are seeking only one particular so- 
lution. 

Example 6.1: For the special case n = 3, Equations 6.9 reduce to 

v'iyi+v' 2 y 2 +v' 3 y 3 =0 

v[y[+v 2 y 2 +v' 3 y 3 =0 (6.10) 

vM+vW + vtf^Mx) 



For the case n = 2, Equations 6.9 become 

v \y\+ v' 2 y 2 =0 
v[y¡ + v 2 y 2 = <j>(x) 



(6.11) 



and for the case n = 1 , we obtain the single equation 

v[ yi =<Kx) (6.12) 

Since^jí^), y 2 (x),..., y n (x) are n linearly independent solutions of the same 
equation L(y) = 0, their Wronskian is not zero (Theorem 4.3). This means 
that the system 6.9 has a nonzero determinant and can be solved unique- 
\yforv¡(x),v 2 (x),...y n (x). 

Scope of the Method 

The method of variation of parameters can be applied to all linear differ- 
ential equations. It is therefore more powerful than the method of unde- 
termined coefficients, which is restricted to linear differential equations 
with constant coefficients and particular forms of 0(x). Nonetheless, in 
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those cases where both methods are applicable, the method of undeter- 
mined coefficients is usually the more efficient and, henee, preferable. 

As a practical matter, the integration of v'¡(x) may be impossible to 
perform. In such an event other methods (in particular, numerical tech- 
niques) must be employed. 

Initial-Value Problems 

Initial- valué problems are solved by applying the initial conditions to the 
general solution of the differential equation. It must be emphasized that 
the initial conditions are applied only to the general solution and not to 
the homogeneous solution y h that possesses all the arbitrary constants that 
must be evaluated. The one exception is when the general solution is the 
homogeneous solution; that is, when the differential equation under con- 
sideration is itself homogeneous. 

Solved Problems 



Solved Problem 6.1 Solve y" -y'-2y = 4x 2 . 

From Problem 5.l f y h = c l e 2K + c 2 e~ x . Here (¡){x) = 4x 2 , a second degree 
polynomial. Using Equation 6.1, we assume that 

y p =A 2 x 2 +A l x + A (6.13) 

Thus, y = 2A 2 x + A, and y" = 2A 2 . Substituting these results into the 
differential equation, we have 

2A 2 - (2A 2 x + A, ) - 2(A 2 x 2 + A, x + A ) = 4x 2 

or, equivalently, 

(-2A 2 )x 2 + (-2A 2 - 2A, )x + (2A 2 - A, - 2A ) = 4x 2 + (0)x + 

Equating the coefficients of like powers of x, we obtain 

-2A 2 =4 -2A 2 -2A,=0 2A 2 -A,-2A =0 
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Solving this system, we find that A 2 = -2, A, = 2, and A Q = -3. Henee 
Equation 6.13 becomes 

y p =-2x 2 +2x-3 

and the general solution is 

y = y h + y p = c 1 e 2x +c 2 e~ x -2x 2 + 2x-3 

Solved Problem 6.2 Solve y" — y' — 2y = sin2x. 

Again by Problem 5.1, y h = c x e 2x +c 2 e~ x . Here (p(x) has the form dis- 
played in Case 3 with k x = 1, k 2 = 0, and ¡3 = 2. Using Equation 6.3, we 
assume that 

y =Aún2x + Bcos2x (6.14) 

Thus, y = 2Acos2x-25sin2xand> ! '' = -4Asin2x-45cos2x. Substi- 
tuting these results into the differential equation, we have 

(-4A sin 2x - 4B eos 2x) - (2 A eos 2x - 2físin 2x) 

-2(A sin 2x + B eos 2x) = sin 2x 

or, equivalently, 

(-6A + 2fi)sin2x + (-6B- 2A)cos2x = (l)sin2x + (0)cos2a; 

Equating coefficients of like terms, we obtain 

-6A + 2B = 1 -2A-6B = 

Solving this system, we find that A = -3/20 and B = 1/20. Then from 
Equation 6.14, 



20 20 



y =-— sin2.\; + — cos2x 



and the general solution is 
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y = c,e 2x +c~,e x sin2xH coslx 

J x 20 20 



Solved Problem 6.3 Solve y'" + y' = sec*. 
This is a third-order equation with 

y h = c x +c 2 cosx + c 3 sin x 
It folio ws from Equation 6.8 that 

y p = v x +v 2 cosx + v 3 sinx (6.15) 

Here y l = 1, y 2 = eos x, y 3 = sin x, and (p(x) = sec x, so Equation 6. 10 be- 
comes 

v[ (1) + v 2 (eos x) + v' 3 (sin x) = 
v[ (0) + v 2 (- sin x) + V3 (eos x) = 
v[ (0) + v' 2 (- eos x) + V3 (- sin x) = sec x 

Solving this set of equations simultaneously, we obtain v[=secx, 
v 2 = -1, and V3 = -tanx. Thus, 

v i - v \dx = seexí/x = ln I sec x + tan x I 
v 2 = f v^í/x = í— ld!x = -x 

v 3 = ( v' 3 dx = ( - tan xdx = - f dx = ln I eos a: I 

i ' J i ros r 



Substituting these valúes into Equation 6.15, we obtain 

y =lnlsecx + tanx l-xcosA: + (sinx)lnlcosA: I 
The general solution is therefore 

y = yh+y P 

= c x +c 2 cosx + c 3 sin x + \n Isecx + tanx l-xcosx + (sinx)ln Icosx I 
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In This Chapter: 

• Spring Problems 

• Eléctrica! Circuit Problems 

• Buoyancy Problems 

• Classifying Solutions 

• Solved Problems 

Spring Problems 

The simple spring system shown in Figure 7-1 consists of a mass m at- 
tached to the lower end of a spring that is itself suspended vertically from 
a mounting. The system is in its equilibrium position when it is at rest. 
The mass is set in motion by one or more of the following means: dis- 
placing the mass from its equilibrium position, providing it with an ini- 
tial velocity, or subjecting it to an external forcé F(t). 

47 
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Equilibrium position 



x = 



Initial position at t = 

-i. 




Positive jr-direction 

Figure 7-1 

Hooke's law: The re storing forcé F of a spring is equal and opposite to 
thef orces applied to the spring and is proportional to the extensión (con- 
traction) l ofthe spring as a result ofthe applied forcé; that is, F = —kl, 
where k denotes the constant of proportionality, generally called the 
spring constant. 



Example 7.1. A steel ball weighing 128 Ib is suspended from a spring, 
whereupon the spring is stretched 2 ft from its natural length. The applied 
forcé responsible for the 2-ft displacement is the weight of the ball, 128 Ib. 
Thus, F = -128 Ib. Hooke's law then gives -128 = -k(2), or k = 64 lb/ft. 

For convenience, we choose the downward direction as the positive 
direction and take the origin to be the center of gravity of the mass in the 
equilibrium position. We assume that the mass of the spring is negligible 
and can be neglected and that air resistance, when present, is proportion- 
al to the velocity of the mass. Thus, at any time t, there are three forces 
acting on the system: (1) F(t), measured in the positive direction; (2) a 
restoring forcé given by Hooke's law as F s = -kx, k>0; and (3) a forcé 
due to air resistance given by F a = —ax, a>0, where a is the constant of 
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proportionality. Note that the restoring forcé F s always acts in a direction 
that will tend to return the system to the equilibrium position: if the mass 
is below the equilibrium position, then x is positive and —kx is negative; 
whereas if the mass is abo ve the equilibrium position, then x is negative 
and —kx is positive. Also note that because a > the forcé F a due to air 
resistance acts in the opposite direction of the velocity and thus tends to 
retard, or damp, the motion of the mass. 

It now follows from Newton's second law (see Chapter Three) that 
m'x = -kx — ax + F(t) , or 

.. a . k F(t) 

x + —x + — x = ^^- (7.1) 

m m m 

If the system starts at t = with an initial velocity v and from an initial 
position x Q , we also have the initial conditions 

x(0) = x o x(0) = v o (7.2) 

The forcé of gravity does not explicitly appear in 7. 1, but it is pres- 
ent nonetheless. We automatically compensated for this forcé by mea- 
suring distance from the equilibrium position of the spring. If one wish- 
es to exhibit gravity explicitly, then distance must be measured from the 
bottom end of the natural length of the spring. That is, the motion of a vi- 
brating spring can be given by 

.. a . k F(t) 

x + — x + —x = g + —^- 
m m m 

if the origin, x = 0, is the terminal point of the unstretched spring before 
the mass m is attached. 



Electrical Circuit Problems 

The simple electrical circuit shown in Figure 7-2 consists of a resistor R 
in ohms; a capacitor C in farads; an inductor L in henries; and an elec- 
tromotive forcé (emf ) E(t) in volts, usually a battery or a generator, all 
connected in series. The current / flowing through the circuit is measured 
in amperes and the charge q on the capacitor is measured in coulombs. 
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-> — WMV- 



£(/) 



w 

Figure 7-2 

Kirchhoff's loop law: The algebraic sum ofthe voltage drops in a sim- 
ple closed electric circuit is zero. 



It is known that the voltage drops across a resistor, a capacitor, and 
an inductor are respecti vely RI, ( 1 IC)q, and L(dlldi) where q is the charge 
on the capacitor. The voltage drop across an emf is —E(t). Thus, from 
Kirchhoff's loop law, we have 



RI+L— + — q-E(t) = 
dt C 

The relationship between q and / is 

dq di _ d 2 q 

dt dt dt 2 

Substituting these valúes into Equation 7.3, we obtain 



1 



1 



d A q Rdq . , 

dt 2 L dt LC L 



The initial conditions for q are 



dt 



= I(0) = I o 



(7.3) 



(7.4) 



(7.5) 



(7.6) 



To obtain a differential equation for the current, we differentiate 
Equation 7.3 with respect to t and then substitute Equation 7.4 directly 
into the resulting equation to obtain 



dt 2 



Rdl 1 , 1 dE(t) 

+ + / = — 

L dt LC L dt 



(7.7) 
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The first initial condition is 7(0) = I Q . The second initial condition is ob- 
tained from Equation 7.3 by solving for clllclt and then setting t = 0. Thus, 



di 

dt 



= L E (o)-^i Q -L ( 

t=Q ^ ^ Lj\s 



(7.8) 



An expression for the current can be gotten either by solving Equation 
7.7 directly or by solving Equation 7.5 for the charge and then differ- 
entiating that expression. 

Buoyancy Problems 

Consider a body of mass m submerged either partially or totally in a liq- 
uid of weight densityp. Such a body experiences two forces, adownward 
forcé due to gravity and a counter forcé governed by: 

Archimedes' principie: A body in liquid experiences a buoyant upward 
forcé equal to the weight ofthe liquid displaced by that body. 



Positive x direction 

4 



Equilibrium siate 



.v = 



Equilibrium position 



T 




- Water line 



Figure 7-3 



Equilibrium occurs when the buoyant forcé of the displaced liquid 
equals the forcé of gravity on the body. Figure 7-3 depicts the situation 
for a cylinder of radius r and height H where h units of cylinder height 
are submerged at equilibrium. At equilibrium, the volume of water dis- 
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placed by the cylinder is nr 2 h, which provides a buoyant forcé of nr 2 hp 
that must equal the weight of the cylinder mg. Thus, 

nr 2 hp = mg (7.9) 

Motion will occur when the cylinder is displaced from its equilibri- 
ura position. We arbitrarily take the upward direction to be the positive 
x-direction. If the cylinder is raised out of the water by x(t) units, as shown 
in Figure 7-3, then it is no longer in equilibrium. The downward or neg- 
ative forcé on such a body remains mg but the buoyant or positive forcé 
is reduced to nr 2 \h - x(t)]p. It now follows from Newton's second law 
that 

mx = nr 2 [h - x(t)]p - mg 

Substituting 7.9 into this last equation, we can simplify it to 

mx = —nr 2 x{t)p 



or 



x+EL± x = o (v.io) 

m 



Classifying Solutions 

Vibrating springs, simple electrical circuits, and floating bodies are all 
governed by second-order linear differential equations with constant co- 
efficients of the form 

x + a x x + a x = f(t) (7.11) 

For vibrating spring problems defined by Equation 7.1, a { = a/m, a Q = 
k/m, and/(í) = F(t)/m. For buoyancy problems defined by Equation 7. 10, 
a 1 = 0, a Q = nr 2 plm, and/(í) = 0. For electrical circuit problems, the in- 
dependent variable x is replaced either by q in Equation 7.5 or / in Equa- 
tion 7.7. 

The motion or current in all of these systems is classified asfree and 
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undamped when/(í) = and a l = 0. It is classi- 

fied asfree and damped when/(í) is identically 

zero but a x is not zero. For damped motion, 

there are three sepárate cases to consider, ac- 

cording as the roots of the associated character- 

istic equation (see Chapter Fi ve) are ( 1 ) real and 

distinct, (2) equal, or (3) complex conjúgate. 

These cases are respectively classified as (1) 

overdamped, (2) critically damped, and (3) os- 

cillatory damped (or, in electrical problems, underdamped). lf/(í) is not 

identically zero, the motion or current is classified as forcé d. 

A motion or current is transient if it "dies out" (that is, goes to zero) 
as t — > oo. A steady-state motion or current is one that is not transient and 
does not become unbounded. Free damped systems always yield transient 
motions, while forced damped systems (assuming the external forcé to be 
sinusoidal) yield both transient and steady-state motions. 

Free undamped motion defined by Equation 7.11 with a { = and 
/(?) = always has solutions of the form 




x(t) = c, eos cot + c 2 sin cot 



(7.12) 



which defines simple harmonio motion. Here c v c 2 , and co are constants 
with co often referred to as circular f re que ncy. The natural frequencyfis 



/=- 

and it represents the number of complete oscillations per time unit un- 
dertaken by the solution. The period of the system of the time required to 
complete one oscillation is 

f 



Equation 7.12 has the altérnate form 

x(t) = (-l) k A cos(cot - (¡)) 



(7.13) 



where the amplitude A = -yc, +c 2 , the p has e angle (j> = arctan(c 2 /c 1 ), 
and k is zero when c. is positive and unity when c. is negative. 
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Solved Problems 

Solved Problem 7.1 A 10-kg mass is attached to a spring, stretching it 
0.7 m from its natural length. The mass is started in motion from the 
equilibrium position with an initial velocity of 1 m/sec in the upward di- 
rection. Find the subsequent motion, if the forcé due to air resistance is 
-90i N. 

Taking g = 9.8m/sec 2 , we have w = mg = 98 N and k = wll = 140 N/m. 
Furthermore, a = 90 and F{t) = (there is no external forcé). Equation 7. 1 
becomes 

x + 9i + 14*; = (7.14) 

The roots of the associated characteristic equation are X x = -2 and X 2 = 
-7, which are real and distinct; henee this problem is an example of over- 
damped motion. The solution of 7.14 is 

x = c l e~ 2t +c 2 e~ 7 ' 

The initial conditions are x(0) = (the mass starts at the equilibrium 
position) and x(0) = -l (the initial velocity is in the negative direction). 
Applying these conditions, we find that c x =-c 2 =-j , so that x = 

j(e~ 7t -e~ 2 '). Note that x — > as t — > °°; thus, the motion is transient. 

Solved Problem 7.2 Determine whether a cylinder of radius 4 in, height 
10 in, and weight 15 Ib can float in a deep pool of water of weight densi- 
ty 62.5 lb/ft 3 . 

Let h denote the length (in feet) of the submerged portion of the 
cylinder at equilibrium. With r = j ft, it follows from Equation 7.9 that 

h = -^L = !^ « 0.688fí = 8.25in 

xr 2 p 7r(^) 2 62.5 

Thus, the cylinder will float with 10 - 8.25 = 1.75 in of length above the 
water at equilibrium. 



r ~i 
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Definition of the Laplace Transform 

Let/(x) be defined for < x < °° and let s denote an arbitrary real vari- 
able. The Laplace transform off(x), designated by either i£ {f(x) } or F(s), 
is 

£{f(x)} = F(s) = ]e- sx f(x)dx (8.1) 

o 

for all valúes of s for which the improper integral converges. Conver- 
gence occurs when the limit 



R 

1 

R^> 



lim ¡e- sx f(x)dx (8.2) 







exists. If this limit does not exist, the improper integral diverges andf(x) 
has no Laplace transform. When evaluating the integral in Equation 8. 1 , 
the variable s is treated as a constant because the integration is with re- 
spect to x. 

The Laplace transforms for a number of elementary functions are 
given in Appendix A. 

Properties of Laplace Transforms 

Property 8.1 (Linearity). If £{f(x)} = F(s) and £{g(x)} = G(s), then 
for any two constants c l and c 2 



2{cJ(x) + c 2 g(x)} = c^{f(x)} + c 2 £{g(x)} (8.3) 

= c l F(s) + c 2 G(s) 

Property 8.2. If !£{f(x)} = F(s), then for any constant a 

&{e ax f{x)} = F{s-a) (8.4) 

Property 8.3. If !£{f(x)} = F(s), then for any positive integer n 

£{x n f(x)} = (-iy-^[F(s)] (8.5) 

ds 
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f(x) 
Propertv 8.4. If ¿6 {/(*)} = F(s) and if lim^-— - exists, then 

x->0 X 
x>0 

2¡^f(x)\ = JF(t)dt (8.6) 

s 

Property 8.5. If £{f(x)} = F(s), then 

2 \]f(t)dt \ = ]¡F(s) (8.7) 

Property 8.6. If/(x) is periodic with period co, that is,/(x + co) =f(x), then 

O) 

¡e- sx f(x)dx 

mW}= ° x _ e . m (8-8) 

Functions of Other Independent Variables 

For consistency only, the definí tion of the Laplace transform and its prop- 
erties, Equations 8.1 through 8.8, are presented for functions of x. They 
are equally applicable for functions of any independent variable and are 
generated by replacing the variable x in the abo ve equations by any vari- 
able of interest. In particular, the counter part of Equation 8. 1 for the 
Laplace transform of a function of t is 

SE{f(t)} = F(s) = ]e- sí f(t)dt 



Definition of the Inverse Laplace Transform 

An inverse Laplace transform of F(s) designated by i£ -1 {F(s)}, is an- 
other function /(x) having the property that !£{f(x)} = F(s). 

The simplest technique for identifying inverse Laplace transforms is 
to recognize them, either from memory or from a table such as in the Ap- 
pendix. If F(s) is not in a recognizable form, then occasionally it can be 
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transformed into such a form by algebraic manipulation. Observe from 
the Appendix that almost all Laplace transforms are quotients. The rec- 
ommended procedure is to first convert the denominator to a form that 
appears in the Appendix and then the numerator. 

Manipulating Denominators 

The method of completing the square converts a quadratic polynomial 
into the sum of squares, a form that appears in many of the denominators 
in the Appendix. In particular, for the quadratic 



as 2 +bs + c =a\s 2 -\ — s\ + c 



= a 


2 b ( 

s +—s + \ 

a y 


-í 

2a J 


+ 


4a 


= a 


Hh 


í b 
c — 

l 4 


2^ 





= a(s + k) 2 + h 2 



where k = bl2a and h = -\jc-(b 1 4a) . 

The method of partial fractions transforms a function of the form 
a(s)/b(s), where both a(s) and b(s) are polynomials in s, into the sum of 
other fractions such that the denominator of each new fraction is either a 
first-degree or a quadratic polynomial raised to some power. The method 
requires only that the degree of a(s) be less than the degree of b(s) (if this 
is not the case, first perform long división, and consider the remainder 
term) and b(s) be factored into the product of distinct linear and quadrat- 
ic polynomials raised to various powers. 

The method is carried out as follows. To each factor of b(s) of the 
form (s - á) m , assign a sum of m fractions, of the form 



A, 



A, 



s-a (s-a) 



- + ••• + - 



A„ 



(s-a)" 



To each factor of b(s) of the form (s 2 + bs + c) p , assign a sum of p frac- 
tions, of the form 



B,s + C, 



B 2 s + C 2 



s 2 +bs + c (s 2 +bs + c) 2 



+ ••• + - 



(s 2 +bs + c) p 
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Here A ., B., and C. (i = 1,2,..., m; /', & = 1,2,..., p) are constants which still 

/'y Ac v ' ' ' ' •/' ' ' ' ' ' 

must be determined. 

Set the original fraction a(s)lb{s) equal to the sum of the new frac- 
tions just constructed. Clear the resulting equation of fractions and then 
equate coefficients of like powers of s, thereby obtaining a set of simul- 
taneous linear equations in the unknown constants A f , B., and C k . Final- 
ly, solve these equations for A v B., and C k . 

Manipulating Numerators 

A factor s - a in the numerator may be written in terms of the factor 5 - 
b, where both a and b are constants, through the identity s - a = (s — b) + 
(b - a). The multiplicative constant a in the numerator may be written ex- 
plicitly in terms of the multiplicative constant b through the identity 

„=£(») 

b 

Both identities genérate recognizable inverse Laplace transforms when 
they are combined with: 

Property 8.7 (Linearity). If the inverse Laplace transforms of two func- 
tions F(s) and G(s) exist, then for any constants c { and c 2 , 

£~ 1 {c 1 F(s) + c 2 G(s)} = c J £~ l [F(s)} + c 2 £- 1 {G(s)} 

Convolutions 

The convolution of two functions/(x) and g(x) is 

X 

f(x)*g(x) = jf(t)g(x-t)dt 



(8.9) 



Theorem 8.1. f(x) * g(x) = g(x) *f(x). 



Theorem 8.2. (Convolution Theorem). If ¿£ {/(*)} = F(s) and £{g(x)} 
= G(s), then £{f(x) * g(x)} = ${f(x)} ${g(x)} = F(s)G(s) 
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You Need to Know 



/ 



The inverse Laplace transform of a product is 
computed using a convolution. 

£-\F(s)G(s)} = f(x)*g(x)=g(x)*f(x) (8.10) 



If one of the two convolutions in Equation 8.10 is simpler to calcúlate, 
then that convolution is chosen when determining the inverse Laplace 
transform of a product. 



Unit Step Function 

The unit step function u(x) is defined as 

Í0 x<0 



u(x) = 



1 x>0 



As an immediate consequence of the definition, we have for any num- 
ber c, 



u(x-c) = 



x<c 

1 x>c 



The graph of u(x - c) is gi ven in Figure 8- 1 . 
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Figure 8-1 
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Theorem 8.3. £{u(x-c)} = -e~ cs . 

s 

Translations 

Given a function/(x) defined forx > 0, the function 

x<c 



u(x — c)f(x — c) = 



f(x -c) x>c 



represents a shift, or translation, of the function /(x) by c units in the pos- 
itive x-direction. For example, if f(x) is given graphically by Figure 8-2, 
then u(x - c)f(x - c) is given graphically by Figure 8-3. 



/<*> 
A 




Figure 8-2 



-»-,' 



■íU-rí/U-e) 

4 




■+■.% 



Figure 8-3 
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Theorem 8.4. If F(s) = <£{/(x)}, then 

£{u{x - c)f(x -c)} = e~ cs F(s) 
Conversely, 



1 {e- cs F(s) } = u(x- c)f(x -c) = 



x<c 

f(x -c) x>c 



Solved Problems 

Solved Problem 8.1 Find £{ e™ }. 
Using Equation 8. 1 , we obtain 



F(s)= ¡e- sx e ox dx= lim ¡e ( °- s)x dx 





= lim 



a (.a-s)x 



a-s 



-i.x=R 



= lim 

R->~ 



Jo-s)R 



-1 



a-s 



s-a 



-(fors > a) 



Note that when s < a, the improper integral diverges. (See also entry 7 in 
the Appendix.) 

Solved Problem 8.2 Find £{xe 4x }. 

This problem can be done three ways. 

(a) Using entry 14 of the Appendix with n = 2 and a = 4, we have 
directly that 

SE{xe 4x } = — L-j. 

(s-4) 2 

(b) Set/(x) = x. Using Property 8.2 with a = 4 and entry 2 of the Ap- 
pendix, we have 
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F(s) = £{f(x)} = ${x} = \ 
s 



and 

£{e 4x x} = F(s-4) = l -^r 

(5-4) 2 

(c) Set/(x) = e 4x . Using Property 8.3 with n = 1 and the results of 
Problem 8. 1, or alternatively, entry 7 of the Appendix with a = 4 
we find that 

4.v i ' 



F(s) = £{f(x)} = ${e 4x } = 



5-4 

and 



4x, -■ Cl ( \ 



£{xe« x } = -F'(s) = 



dsys-4J (s-4y 

s + 3 



Solved Problem 8.3 Use partial fractions to decompose 

(s-2)(s + \) 

To the linear factors s - 2 and s + 1 , we associate respectively the 
fractions A /(s - 2) and Bl{s + 1). We set 

5 + 3 A B 

- + - 



(s-2)(s + l) s-2 5 + 1 
and, upon clearing fractions, obtain 

s + 3 = A(s + \) + B(s-2) (8-11) 

To find A and B, we substitute s = -1 and 5 = 2 into 8. 1 1 , we immediate- 
ly obtain A = 5/3 and B = -2/3. Thus, 

5 + 3 5/3 2/3 



Solved Problem 8.4 Find 



(5-2)(5 + l) 5-2 5 + 1 

-I I 5 + 3 



(5-2)(5 + l) 
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No function of this form appears in the Appendix. Using the results 
of Problem 8.3 and Property 8.7, we obtain 



5 + 3 1 5.._, í 1 1 2„_, f 1 



(s-2)(s + l)l 3 ls-2) 3 U + l 



5 ix 2 x 
= —e — e 
3 3 
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In This Chapter: 

• Laplace Transforms of Derivatives 

• Solutions of Linear Differential 
Equations with Constant 
Coefficients 

• Solutions of Linear Systems 

• Solved Problems 

Laplace Transforms of Derivatives 

Denote !£{y(x)} by Y(s). Then under very broad conditions, the Laplace 
transform of the «th-derivative (n = 1,2,3»—) ofy(x) is 



d n y\ = s"Y(s)-s"- l y(0)-s n - 2 y'(0)- 
dx n \ ■■■-sy (n - 2 \0)-y (n - l \0) 



(9.1) 



If the initial conditions on y(x) at x = are given by 

y(0) = c o ,y'(0) = c l ,...,y ( "- i) (0) = c n _ l (9.2) 
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then (9.1) can be rewritten as 

gl = s"Y( S )-c s"-' - c,5"- 2 — - c„_ 2 s - c„_, (9.3) 
For the special cases of n = 1 and n = 2, Equation 9.3 simplifies to 

<Z{yix)) = sY(s)-c (9.4) 



£{y"(x)} = s 2 Y(s)-c s- Cl 



(9.5) 



ik 



Note! 



Laplace transforms convert differential equations 
into algebraic equations. 



Solutions of Linear Differential Equations 
with Constant Coefficients 

Laplace transforms are used to solve initial- valué problems given by the 
nth-order linear differential equation with constant coefficients 



fc d "y h 



d n ~ x y 

dx n ~ x 



. +b y +by=g ( X ) 

dx 



(9.6) 



together with the initial conditions specified in Equation 9.2. First, take 
the Laplace transform of both sides of Equation 9.6, thereby obtaining an 
algebraic equation for Y(s). Then solve for Y(s) algebraic ally, and final- 
ly take inverse Laplace transforms to obtain y{x) = i£ -1 { Y(s)}. 

Unlike previous methods, where first the 
differential equation is solved and then the 
initial conditions are applied to evalúate 
the arbitrary constants, the Laplace trans- 
form method sol ves the entire initial- valué 
problem in one step. There are two excep- 
tions: when no initial conditions are speci- 
fied and when the initial conditions are not 
at x = 0. In these situations, c Q through c nl 
in Equations 9.2 and 9.3 remain arbitrary 
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and the solution to differential equation 9.6 is found in terms of these con- 
stants. They are then evaluated separately when appropriate subsidiary 
conditions are provided. 

Solutions of Linear Systems 

Laplace transforms are useful for solving systems of linear differential 
equations; that is, sets of two or more differential equations with an equal 
number of unknown functions. If all of the coefficients are constants, then 
the method of solution is a straightforward generalization of the one de- 
scribed above. Laplace transforms are taken of each differential equation 
in the system; the transforms of the unknown functions are determined 
algebraically from the resulting set of simultaneous equations; inverse 
transforms for the unknown functions are calculated with the help of the 
Appendix. 

Solved Problems 

Solved Problem 9.1 Solve / - 5v = e 5x ; v(0) = 0. 

Taking the Laplace transform of both sides of this differential equation 
and using Property 8. 1 , we find that #{/} - 5i£ {y} = £{e 5x }. Then, us- 
ing the Appendix and Equation 9.4 with c Q = 0, we obtain 

[sY(s) - 0] - 5Y(s) = from which Y(s) = - 

s-5 (-^-5) 

Finally, taking the inverse Laplace transform of Y(s), we obtain 
y(x) = i£-' {Y(s) } = %- 1 j^-rl = xe 5x 



(s-5Y 



(see Appendix, entry 14). 

Solved Problem 9.2 Solve the system 

y" + z + y = 



y(0) = 0, y'(0) = 0, z(0) = \ 



68 DIFFERENTIAL EQUATIONS 

Denote !£{y(x)} and !£{z(x)} by Y(s) and Z(s) respectively. Then, 
taking the Laplace transforms of both differential equations, we obtain 



[s 2 Y(s) - (0)s - (0)] + Z(s) + Y(s) = 
[sZ(s)-l] + [sY(s)-0] = 



or 



(s 2 +l)Y(s) + Z(s) = 

Y(s) + Z(s) = - 

s 



Solving this last system for Y(s) and Z(s), we find that 

Y(s) = -\ Z(s) = - + \ 
s s s 

Thus, taking inverse transforms, we conclude that 
y(x) = --x 2 z(x) = \ + -x 2 
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Matrices and Vectors 

A matrix (designated by an uppercase boldface letter) is a rectangular ar- 
ray of elements arranged in horizontal rows and vertical columns. In this 
book, the elements of matrices will always be numbers or functions of 
the variable t. If all the elements are numbers, then the matrix is called a 
constant matrix. 

Example 10.1. 



[1 2] 


1 é 2 


, and 


3 4j 


t -1 1 





1 t l 



cosí 



are all matrices. In particular, the first matrix is a constant matrix, where- 
as the last two are not. 

A general matrix A having p rows and n columns is given by 









a u 


a ]2 ■ 


■ <hn 




r i 




«21 


a 22 ■ 


■ a 2n 


A = 


H 
















a,A 


a p2 ■ 


■ a pn 



where a . . represents that element appearing in the /th row andy'th column. 



1^ 



Note! 



A matrix is square if it has the same number of rows 
and columns. 



A vector (designated by a lowercase boldface letter) is a matrix having 
only one column or one row. (The third matrix given in Example 10. 1 is 
a vector.) 
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Matrix Addition 

The sum A + B of two matrices A = [a¡ .] and B = [b¡ ] having the same 
number of rows and the same number of columns is the matrix obtained 
by adding the corresponding elements of A and B. That is, 

A + B = [a ¡J ] + [b u ] = [a u +b u ] 

Matrix addition is both associative and commutative. Thus, A + (B + C) 
= (A + B) + C and A + B = B + A. 

Scalar and Matrix Multiplication 

If X is either a number or a function of t, then XA (or, equivalently, Al) is 
defined to be the matrix obtained by multiplying every element of A by 
X. That is, 

XA = X[a¡j] = [Xa ij ] 

Let A = [a.] and B = [b¡] be two matrices such that A has r rows and n 
columns and B has n rows and p columns. Then the product AB is de- 
fined to be the matrix C = [c ..] given by 



c ü 



= Y. a ¡k b kj 0' = U,. ..,/-;; = 1,2,. ..,/?) 



k=\ 



The element c ¡ . is obtained by multiplying the elements of the z'th row of 
A with the corresponding elements of they'th column of B and summing 
the results. 

Matrix multiplication is associative and distributes over addition; 
in general, however, it is not commutative. Thus, 

A(BC) = (AB)C, A(B + C) = AB + AC, and (B + C)A = BA + CA. 



Important! 

In general, AB * BA. 




72 DIFFERENTIAL EQUATIONS 

Powers of a Square Matrix 

If n is a positive integer and A is a square matrix, then 

A W =AA-A 

n times 



In particular, A 2 = A A and A 3 = AAA. By definí tion, A° = I, where 



1 = 



[1 





• 


• 


01 





1 


• 


• 











1 • 


• 











• 


• 1 











• 


• 


1 



is called an identity matrix. For any square matrix A and identity matrix 
I of the same size 

AI = IA = A 



Differentiation and Integration of Matrices 

The derivative of A = [a..] is the matrix obtained by differentiating each 
element of A; that is, 



dX 
dt 



dt 



Similarly, the integral of A, either definite or indefinite, is obtained by in- 
tegrating each element of A. Thus, 



\Adt = 



\a-.dt 



and Adt = 



JV 
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The Characteristic Equation of a Matrix 

The characteristic equation of a square matrix A is the polynomial equa- 
tion in X given by 

det(A-AI) = (10.1) 

where det( ) stands for "the determinant of." Those valúes of X which sat- 
isfy 10.1, that is, the roots of 10.1, are the eigenvalues of A, a^-fold root 
being called an eigenvalue of multiplicity k. 

Theorem 10.1. {Cayley-Hamilton theorem). Any square matrix satis- 
fies its own characteristic equation. That is, if 

det(A - XI) = b n X" + b n _ x X n ~ x +■■■ + b 2 X 2 + b x X + b 

then b„A" + b n _ x A" -1 + • • • + b 2 A 2 + b x A + b I = 0. 

Definition of the Matrix Exponential e AÍ 

For a square matrix A, 

e As =I + -At + -A 2 1 2 +-=J J -A"t H (10.2) 

1! 2! ¿tn\ 

n=0 

The infinite series 10.2 converges for every A and t, so that e At is defined 
for all square matrices. 

Computation of the Matrix Exponential e AÍ 

For actually computing the elements of e At , 10.2 is not generally useful. 
However, it follows (with some effort) from Theorem 10. 1, applied to the 
matrix Ai, that the infinite series can be reduced to a polynomial in t. 
Thus: 

Theorem 10.2. If A is a matrix having n rows and n columns, then 

e At = a n _ x A n - x t n ~ x + a n _ 2 A"- 2 t n ~ 2 +■■■ + a x At + a l (10.3) 
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where or n , a.,...,a _ , are functions of t which must be determined for 



"0' "1'"''"b - 1 



each A. 

Example 10.2. When A has two rows and two columns, then n = 2 and 



e At =or,Aí + or I 



When A has three rows and three columns, then n = 3 and 



e At = a,A 2 t 2 +a,At + or n I 



(10.4) 



(10.5) 



Theorem 10.3. Let A be as in Theorem 10.2, and define 



r(X) = a n _ x X" ] +a n _ 2 X" 2 + --- + a 2 X 2 +a x X + a 



Then if X ¡ is an eigen valué of Ai, 



e x, =T{X t ) 



(10.6) 



(10.7) 



Furthermore, if X¡ is an eigen valué of multiplicity k, k > 1 , then the fol- 
lowing equations are also valid: 



e x >=±r{X) 

dX 



di 2 



.k-\ 

A-, a í 1 \ 

6 = ^ rW 



(10.8) 



Note that Theorem 10.3 involves the eigenvalues of Ai; these are t 
times the eigenvalues of A. When computing the various derivatives in 
10.8, one first calculates the appropriate derivatives of the expression 
10.6 with respect to X, and then substitutes X = X r The reverse procedure 
of first substituting X = X¡ (a function of i) into 10.6, and then calculating 
the derivatives with respect to t can give erroneous results. 
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Example 10.3. Let A have four rows and four columns and let X = 5t and 
X = 2t be eigen valúes of Ai of multiplicities three and one, respectively. 
Then n = 4 and 

r{X) = a 3 X 3 + a 2 X 2 + a x X + a 

r'(X) = 3a 3 X 2 + 2a 2 X + a x 
r"{l) = 6a 3 X + 2a 2 



Since X = 5t is an eigenvalue of multiplicity three, it follows that e 5r = 
r(5t), e 5t = r'(5t) and e 5t = r"(5t). Thus, 

e 51 = a 3 (5tf + a 2 (5tf + a { (5?) + a 
e 5 ' =3a 3 (5t) 2 +2a 2 (5t) + a ] 
e 5 ' =6a 3 {5t) + 2a-, 



Also, since X = 2t is an eigenvalue of multiplicity one, it follows that e 2 ' 
= r(2t), or 

e 2 ' =a 3 (2í) 3 + a 2 (2t) 2 + a x (2t) + a 

Notice that we now have four equations in four unknown a's. 

Method of computation: For each eigenvalue X¡ of A/, apply Theorem 
10.3 to obtain a set of linear equations. When this is done for each eigen- 
value, the set of all equations so obtained can be solved for a ,a i ,..., 
a n _ r These valúes are then substituted into Equation 10.3 which, in turn, 
is used to compute e At . 

Solved Problems 

Solved Problem 10.1 Find AB and BA for 



A = 



1 2 3 
4 5 6 



B = 



7 0' 

8 -1 



Since A has three columns and B has two rows, the product AB is not de- 
fined. But 
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BA = 


7 

8 -1 


1 
4 


2 
5 


3 

6 










r7(l) + (0)(4) 


7(2) +(0X5) 7(3)+ (0X6)1 




8(l) + (-l)(4) 


8(2) + (-1)(5) 8(3) + (-l)(6)_ 




[7 14 21] 








4 11 18 






Solved Problem 10.2 Find the 


eigenvalues of A = 


[1 31 
4 2 




We have 

















A-AI = 



+ (-A) 



+ 



Henee, 



det(A - AI) = det 



1 
1_, 

-A l_p-A 3 
-A _~ 4 2-A 
1-A 3 

4 2-A_ 

= (1 - A)(2 - A) - (3)(4) = A 2 - 3A - 10 



The characteristic equation of A is A 2 - 3A -10 = 0, which can be factored 
into (A - 5)(A + 2) = 0. The roots of the equation are Aj = 5 and A 2 = -2, 
which are the eigenvalues of A. 

"1 ll 



Solved Problem 10.3 Find e At for A = 
Here n = 2. From Equation 10.4, 



9 1 



.Ai 



= a,At + a ñ l = 



a,í + a, 

9a,í 



o 






(10.9) 



and from Equation 10.6, r(X) = g^A + a Q . The eigenvalues of Ai are A x = 
4í and A 2 = -2í, which are both of multiplicity one. Substituting these val- 
úes successively into Equation 10.7, we obtain the two equations 



e 4t = Ata, + a n 



e 2t = -2ta, + a n 



CHAPTER 10: Matrices and the Matrix Exponential 77 

Solving these equations for a x and a , we find that 

a x =-{e At -e- 2t ) and a =-(e 4t +2e~ 2t ) 
6í 3 

Substituting these valúes into 10.9 and simplifying, we have 



e 4 '- 1 


[3e 4t +3e 


-2í 


e 4t -i 


,-2' 






e — — 

6 


9e 4t -9e~ 2t 3e 4t +3e~ 2t \ 




Solved Problem 10.4 Find é? A? for A = 


ro 1 1 

8 -2 




Since n = 2, it folio ws from Equations 10.4 and 10.6 that 


e Al =a,Aí + a I = 


[ a a x t 
%a x t -2a { t + a 





(10.10) 



and r{X) = a x X + «q. The eigen valúes of Ai are Á l = 2t and A 2 = -4t, which 
are both of multiplicity one. Substituting these valúes successively into 
Equation 10.7, we obtain 

e 2t =a l (2t) + a 
e~ 4 ' = a x (-4í) + a 

Solving these equations for a x and a , we find that 

a { =— (e 2t -e- 4 ') and a =-(2e 2t +e~ 4 ') 
6t 3 

Substituting these valúes into 10.10 and simplifying, we have 



6 



4£r'+2í>- 



2^ 2 '+4e~ 
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Reduction of Linear Differential Equations 
to a First-Order System 

Reduction of One Equation 

Every initial- valué problem of the form 

b nW^ + b n _ í (t)^ + .~ + b 1 (t)x + b (t) X = g(t>, (11.1) 

x(t ) = c , x(t ) = c l ,...,x (n - l) (t ) = c n _ l (11.2) 

with b n (t) i¿ 0, can be reduced to the first-order matrix system 

x(t) = A(t)x(t) + f(t) 
x(t ) = c 

where A(í), f(t), c, and the initial time t Q are known. The method of re- 
duction is as follows. 

Step 1. Rewrite 11.1 so that d n x/df appears by itself. Thus, 

d n x _ d n ~ l x 

df " a "- l(í) dt"- 1 (11.4) 

+ --- + a í (t)x + a (t)x + f(t) 

where aft) = -bft)/b n (t) (j = 0,1,..., n - 1) and/(í) = g(t)/b n (t). 

Step 2. Define n new variables (the same number as the order of the orig- 
inal differential equation), x^t), x 2 (t),..., x n (t), by the equations 

x l (t) = x(t), x 2 (t) = ^l, 

d < (11.5) 

x(t)= ^m x(t) - d "~ ,x ^ 
3(0 dt 2 -■' n(> df-' 

These new variables are interrelated by the equations 
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x 1 (t) = x 2 (t) 
x 2 (t) = x 3 (t) 
x 3 (t) = x 4 (t) 



(11.6) 



Vi (0 = *„(0 



Step 3. Express cNdt in terms of the new variables. Proceed by first dif- 
ferentiating the last equation of 11.5 to obtain 



. , , d 

x " {t) = T* 
dt 



di"'* 



d n x(t) 
dt" 



Then, from Equations 11.4 and 11.5, 

x n (t) = a n _ x (t) ¿p + • • • + flj (t)x(t) + a (t)x(t) + f(t) 
= a n _ x (t)x n (t) + --- + a l (t)x 2 (t) + a (t) X] (t) + f(t) 

For convenience, we rewrite this last equation so thatXj(í) appears before 
x 2 (t), etc. Thus, 

x n (t) = a (t) X] (?) + o, (t)x 2 (0 + • • • + a n -\ (0*„ (t) + f(t) (11.7) 

Step 4. Equations 11.6 and 11.7 are a system of first-order linear differ- 
ential equations in x x (t), x 2 (t),..., xjt). This system is equivalent to the 
single matrix equation \(t) = A(t)x(f) + f(í) if we define 



*(0 = 



x 2 (t) 



(11.8) 



f(í) = 






o 

f(t)\ 



(11.9) 
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A(í) = 






1 





•• 











1 


•• 














1 








a (t) a x {f) a 2 {t) a 3 (t) 



fl»-i(0 



(11.10) 



Step 5. Define 



c = 



Then the initial conditions (11.2) can be given by the matrix (vector) 
equation x(? ) = c. This last equation is an immediate consequence of 
Equations 11.8, 11.5, and 11.2, since 



x(r ) = 



U('o)] 




x(t ) 




\ c ° 1 


* 2 ('o) 


= 


x(t ) 


= 


c \ 


_*»(*o)_ 




x { - l) (t )_ 




_ C n-\ _ 



= c 



Observe that if no initial conditions are prescribed, Steps 1 through 4 by 
themselves reduce any linear differential equation 11.1 to the matrix 
equation x(r) = A(t)x(t) + f(í). 

Reduction of a System 

A set of linear differential equations with initial conditions also can be re- 
duced to system 1 1.3. The procedure is nearly identical to the method for 
reducing a single equation to matrix form; only Step 2 changes. With a 
system of equations, Step 2 is generalized so that new variables are de- 
fined for each of the unknown functions in the set. 
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Solution of the Initial-Value Problem 

By the procedure described abo ve, any initial- valué problem in which the 
differential equations are all linear with constant coefficients, can be re- 
duced to the matrix system 

x(í) = Ax(í) + f(í); x(í ) = c (11.11) 

where A is the matrix oíconstants. The solution to Equation 11.11 is 

x(í) = e M '-'° } c + e At j e~ As f(s)ds (11.12) 



or equivalently 



s (t) = e M '~'° »c + j e M '- s) f (s)ds (11.13) 



In particular, if the initial- valué problem is homogeneous [i.e., f(t) = 0], 
then both equations 11.12 and 11.13 reduce to 

x(t) = e M '~ to) c (11-14) 

In the abo ve solutions, the matrices e (t ~ to \ e~ As , and e A(J ~ s) are easily 
computed from e At by replacing the variable t by t - t Q , -s, and t - s, re- 
spectively. Usually x(t) is obtained quicker from 11.13 than from 11.12, 
since the former equation involves one less matrix multiplication. How- 
ever, the integráis arising in 11.13 are generally more difficult to evalú- 
ate than those in 11.12. 

Solution with No Initial Conditions 

If no initial conditions are prescribed, the solution of x(t) = Ax(í) + f(í) is 
x(t) = e A 'k + e A, je- At f(t)dt (11.15) 

or, when f(í) = 0, 
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x(t) = e At k (11.16) 

where k is an arbitrary constant vector. All constants of integration can 
be disregarded when computing the integral in Equation 11.15, since they 
are already included in k. 

Solved Problems 

Solved Problem 11.1 Put the initial- valué problem 

x + 2x-8x = 0; x(l) = 2, ¿(1) = 3 

into the form of System 1 1 .3. 

Following Step 1 , we write x = —2x + Sx; henee a i (t) = -2, a Q {i) = 8, and 
f(f) = 0. Then, defining x x (t) = x and x 2 {t) = x (the differential equation is 
second-order, so we need two new variables), we obtain x { = x T Follow- 
ing Step 3, we find 

d 2 x 
x-, = — r- = —2x + Sx = -2x + 8x, 
2 dt 2 2 l 



Thus, 



i, = 0x x + lx 2 



x 2 = 8x, - 2x 2 

These equations are equivalent to the matrix equation x(í) = A(t)x(t) + í(t) 
if we define 



x(0 = 



x 2 {t) 



A(í) = 



1 " 

8 -2 



f(0 = 



The differential equation is then equivalent to the matrix equation x(t) = 
A(í)x(í) + f(í), or simply x(í) = A(í)x(í), since f(í) = 0. The initial condi- 

r 2 ~i 

tions can be given by x(í ) = c, if we define t Q = 1 and c = 
Solved Problem 11.2 Solve Je + 2x - &c = 0; x{\) = 2, x(l) = 3. 
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From Problem 11.1, this initial- valué problem is equivalent to Equation 
11.11 with 



x(0 = 



Xi (t) 

x 2 (t) 



A(í) = 



1 " 

8 -2 



f (t) = c = 



í =l 



The solution to this system is given by Equation 11.14. For this A, e At is 
given in Problem 10.4; henee, 



e Mt-t ) =é ,A(f-l) _ 



4e 2(, - l) + 2e- 4(, - l) e 2(, ~ X) - e~ 4 ^ 



oW-V -S^C" 1 ) 



le 



2(?-l) 



+ 4e 



-4(/-l) 



Therefore, 

x(í) = ^ A(í " 1) 
1 



4e 2(, - l) +2e- 4(, - l) e 2{ " l) - e' 4 ^ 



a 2('-l>_Q^('-l> 



2e 



2(f-l) 



+ 4¿ 



-4(/-l) 



' 2(4e 2(t ~ i) + 2e- 4( '- l) ) + 3(e 2(t - l) -e- 4(t - i} ) ' 
2(Se 2( " l) -Se- 4{ '- l) ) + 3(2e 2( '- i) +4e- 4( '- l) ) 



6 6 

22 



and the solution to the original initial- valué problem is 



11 



1 



(f) = x.(f) = — e* t - v >+-e« t - v > 
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Second-Order Linear Equations 
with Variable Coefficients 

A second-order linear differential equation 

b 2 (x)y" + b l (x)y' + b (x)y = g(x) (12.1) 

has variable coefficients when b 2 (x), b x {x), and b Q (x) are not all constants 
or constant múltiples of one another. If b 2 (x) is not zero in a given inter- 
val, then we can divide by it and rewrite Equation 12.1 as 

y" + P(x)y' + Q(x)y = <p(x) (12.2) 

where P(x) = b x {x)lb 2 {x), Q(x) = b Q (x)/b 2 (x), and <j)(x) = g(x)/b 2 (x). In this 
chapter, we describe procedures for solving many equations in the form 
of 12.1 and 12.2. These procedures can be generalized in a straightfor- 
ward manner to solve higher-order linear differential equations with vari- 
able coefficients. 

Analytic Functions and Ordinary Points 

A function/(.x:) is analytic at x Q if its Taylor series about x Q , 

f f { "\x )(x-x r 

converges tof(x) in some neighborhood of x Q . 





You Need to Know! 


Polynomials, sin x, eos x, and e x are analytic every- 


where. 





Sums, differences, and producís of polynomials, sin x, eos x, and e* are 
also analytic every where. Quotients of any two of these functions are an- 
alytic at all points where the denominator is not zero. 
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The point x Q is an orclinary point of the differential equation 12.2 if 
both P(x) and Q(x) are analy tic at x Q . If either of these functions is not an- 
alytic at x Q , then x Q is a singular point of 12.2. 

Solutions Around the Origin of Homogeneous 
Equations 

Equation 12.1 is homogeneous when g(x) = 0, in which case Equation 
12.2 specializes to 

y" + P(x)y' + Q(x)y = (12.3) 

Theorem 12.1. If x = is an ordinary point of Equation 12.3, then the 
general solution in an interval containing this point has the form 



y = ^a n x" =a y l (x) + a x y 2 (x) (12.4) 

7Í=0 

where a Q and a { are arbitrary constants and y x (x) and y 2 (x) are linearly in- 
dependent functions analytic at x = 0. 

To evalúate the coefficients a in the solution furnished by Theorem 12. 1 , 
use the following fi ve-step procedure known as the power series method. 

Step 1. Substitute into the left side of the homogeneous differential 
equation the power series 



y = V a n x n = a + a x x + a 2 x 2 + a 3 x 3 + a 4 x 4 + ■ 



(12.5) 



+ax n +a^x" +l +a.^x n+2 +- 



together with the power series for 

y' = a, +2a 2 x + 3a 3 x 2 +4a 4 x 3 H — 

+na n x n - 1 +(n + \)a n+i x"+(n + 2)a n+2 x n+i +• 

and 



(12.6) 
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y" = 2a 2 + 6a 3 x + \2a 4 x 2 H — 

+n(n-l)a n x n ~ 2 +(n + l)(n)a n+1 x K 

+(n + 2)0 + l)a n+2 x n +■■■ 



(12.7) 



Step 2. Collect powers of x and set the coefficients of each power of 
x equal to zero. 

Step 3. The equation obtained by setting the coefficient of x' 1 to zero 
in Step 2 will contain a . terms for a finite number of y valúes. Solve 
this equation for the a . term having the largest subscript. The re- 
sulting equation is known as the recurrence formula for the given 
differential equation. 

Step 4. Use the recurrence formula to sequentially determine a (j = 
2,3,4,...) in terms of a Q and a v 

Step 5. Substitute the coefficients determined in Step 4 into Equa- 
tion 12.5 and rewrite the solution in the form of Equation 12.4. 

The power series method is only applica- 
ble when x = O is an ordinary point. Although 
the differential equation must be in the form 
of Equation 12.2 to determine whether x = O is 
an ordinary point, once this condition is veri- 
fied, the power series method can be used on 
eitherform 12.1 or 12.2. If P(x)or Q(x)in 12.2 
are quotients of polynomials, it is often sim- 
pler first to multiply through by the lowest 
common denominator, thereby clearing frac- 
tions, and then to apply the power series method to the resulting equation 
in the form of Equation 12.1. 




Solutions Around the Origin of 
Nonhomogeneous Equations 

If (p(x) in Equation 12.2 is analytic at x = O, it has a Taylor series expan- 
sión around that point and the power series method given above can be 
modified to solve either Equation 12.1 or 12.2. In Step 1, Equations 12.5 
through 12.7 are substituted into the left si de of the nonhomogeneous 
equation; the right side is written as a Taylor series around the origin. 
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Steps 2 and 3 change so that the coefficients of each power of x on the 
left side of the equation resulting from Step 1 are set equal to their coun- 
terparts on the right side of that equation. The form of the solution in Step 
5 becomes 

y = a y 1 (x) + a l y 2 (x) + y 3 (x) 

which has the form specified in Theorem 4.4. The first two terms com- 
prise the general solution to the associated homogeneous differential 
equation while the last function is a particular solution to the nonhomo- 
geneous equation. 

Initial-Value Problems 



Vy Important! 

Solutions to initial-value problems are obtained by 
first solving the given differential equation and then 
applying the specified initial conditions. 



Solutions Around Other Points 

When solutions are required around the ordinary point x Q t- 0, it general- 
ly simplifies the algebra if x Q is translated to the origin by the change of 
variables t = x - x Q . The solution of the new differential equation that re- 
sults can be obtained by the power series method about t = 0. Then the 
solution of the original equation is easily gotten by back-substitution. 

Regular Singular Points 

The point x Q is a regular singular point of the second-order homogeneous 
linear differential equation 

y" + P(x)y' + Q(x)y = (12.8) 
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if x is not an ordinary point but both (x - x Q )P(x) and (x - x ) 2 Q(x) are 
analytic atx Q . We only consider regular singular points atx = 0; if this is 
not the case, then the change of variables t = x - x Q will transíate x Q to the 
origin. 

Method of Frobenius 

Theorem 12.2. If x Q = is a regular singular point of 12.8, then the equa- 
tion has at least one solution of the form 

y = x x f j a„x" 

n=0 

where/l and a (n = 0,1,2,...) are constants. This solution is valid in an in- 

n \ -lili 

terval < x < R for some real number R. 

To evalúate the coefficients a n and X in Theorem 12.2, one proceeds as in 
the power series method described above. The infinite series 



n=0 n=0 

= a x x +a { x Á+l +a 2 x x+2 +• 



(12.9) 



with its derivatives 



Y = Xa x x ~ x +(X + \)a x x x +{X + 2)a 2 x x+i +■■■ 

+(A + n- \)a n _ x x x+n - 2 + (X + n)a n x x+n -' (12. 10) 

+(¿ + n + l)a n+l x x+n +- 



and 



y" = l(X - \)a x x ~ 2 +{X + \){X)a x x x ~ x + (X + 2)(X + \)a 2 x x 

+ --- + (X + n-\)(X + n-2)a n _ x x x+ "- 3 (12.11) 

+(A + n){X + n - l)a n x x+n ~ 2 + (X + n + 1)(A + n)a n+l x x+ "~ l +■■■ 
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are substituted into Equation 12.8. Terms with like powers of x are col- 
lected together and set equal to zero. When this is done for x n the result- 
ing equation is a recurrence formula. A quadratic equation in X, called the 
indicial equation, arises when the coefficient of x° is set to zero and a Q is 
left arbitrary. 

The two roots of the indicial equation can be real or complex. If com- 
plex they will occur in a conjúgate pair and the complex solutions that 
they produce can be combined (by using Euler's relations and the identi- 
ty x a±l = x°e ±l lnx ) to form real solutions. In this book we shall, for 
simplicity, suppose that both roots of the indicial equation are real. Then, 
if X is taken as the larger indicial root, X = X x >X 2 , the method of Frobe- 
nius always yields a solution 

y l (x) = x^f 4 a n (X l )x n (12.12) 

»=o 

to Equation 12.8. [We have written a n {X x ) to indicate the coefficients pro- 
duced by the method when X = X y ] 

If P(x) and Q(x) are quotients of polynomials, it is usually easier first 
to multiply 12.8 by their lowest common denominator and then to apply 
the method of Frobenius to the resulting equation. 

General Solution 

The method of Frobenius always yields one solution to Equation 12.8 of 
the form 12. 12. The general solution (see Theorem 4.2) has the form y = 
c i^ i ( x ) + c ^^. x ) where c { and c 2 are arbitrary constants and y 2 (x) is a sec- 
ond solution of 12.8 that is linearly independent from^C*). The method 
for obtaining this second solution depends on the relationship between 
the two roots of the indicial equation. 

Case 1. If X x - X 2 , is not an integer, then 

y 2 (x) = x^-f j a n (X 2 )x n (12.13) 

where y 2 {x) is obtained in an identical manner d&y x (x) by the method 
of Frobenius, using X 2 in place of X { . 
Case 2. If,^ =A 2 ,then 

y 2 (x) = y l (x)\nx + x^-f j b n (X l )x n (12.14) 



92 DIFFERENTIAL EQUATIONS 

To genérate this solution, keep the recurrence formula in terms of X 
and use it to find the coefficients a (n > 1 ) in terms of both X and a„, 

n v ' O' 

where the coefficient a Q remains arbitrary. Substitute these a n into 
Equation 12.9 to obtain a function y(X, x) which depends on the vari- 
ables X and x. Then 



y¿*)= MX ' x) 



Case 3. If X l - X 2 = N, a positive integer, then 



(12.15) 



y 2 (x) = d_ l y l (x)\nx + x Á2 ^d n (X 2 )x n (12.16) 

«=0 

To genérate this solution, first try the method of Frobenius with X T 
If it yields a second solution, then this solution is y 2 (x), having the 
form of 12.16 with d_ x = 0. Otherwise, proceed as in Case 2 to gen- 
érate y(X, x), whence 

y 2 (x) = ^-[(X-X 2 )y(X,x)]\ x=h (12.17) 



Solved Problems 

Solved Problem 12.1 Determine whether x = is an ordinary point of 
the differential equation 

v"-a-/ + 2}> = 

Here P(x) = -x and Q{x) = 2 are both polynomials; henee they are 
analytic every where. Therefore, every valué of x, in particular x = 0, is 
an ordinary point. 

Solved Problem 12.2 Find a recurrence formula for the power series so- 
lution around x = for the differential equation given in Problem 12.1. 

It folio ws from Problem 12. 1 that x = is an ordinary point of the given 
equation, so Theorem 12.1 holds. Substituting Equations 12.5 through 
12.7 into the left si de of the differential equation, we find 
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[2a 2 + 6a 3 x + \2a 4 x 2 +••• + «(«- X)a n x n ~ 2 +(n + \)(n)a n+l x"~ l +■■■] 
-x[a } + 2a 2 x + 3a 3 x 2 +■■■ + na n x n ~ x + (n + \)a n+l x" +■■■] 
+2[a + a x x + a 2 x 2 + ¿z 3 x 3 + • • • + a n x n + a n+x x n+x +•••] = 

Combining terms that contain like powers of x, we have 

(2a 2 +2a ) + x(6a 3 + a l ) + x 2 (\2a 4 ) + x 3 (20a 5 -a 3 ) 

+ --- + x"[(n + 2)(n + l)a n+2 -na n +2a n ] + --- 

= + 0x + 0x 2 +0x 3 +--- + 0x" +■■■ 

The last equation holds if and only if each coefficient in the left-hand side 
is zero. Thus, 

2a 2 +2a =0, 6a 3 +£7j=0, \2a 4 =0, 20a 5 —a 3 = 0, ... 

In general, 

(n + 2)(n + l)a n+2 -(n-2)a n =0, or, 
(n-2) 
a " +2 (n + 2)(n + l) an 

which is the recurrence formula for this problem. 

Solved Problem 12.3 Find the general solution near x = of y" - xy' + 

2y = 0. 

Successively evaluating the recurrence formula obtained in Problem 12.2 
for n = 0,1,2,..., we calcúlate 

a, = -a n 



a 3 =--a x 

c ' 4 = (12.18) 

a, = — a 3 - — — a, \ = a, 

5 20 3 20 1 6 V 120 ' 
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2 l ,m 

a. = — a, = — (0) 

6 30 4 15 



\ ( 1 } 1 



a 7 = — a* = — a, = a, 

1 42 5 14 1 120 l ) 1680 ' 

4 1 

a„=—a ñ =— (0) = 
8 56 6 14 



Note that since a 4 = 0, it follows from the recurrence formula that all the 
even coefficients beyond a 4 are also zero. Substituting 12.18 into Equa- 
tion 12.5 we have 

v = a n + a.x - a a x 2 — a x x 3 +0x 4 a.x 5 +0x 6 a.x 1 

y ° ' ° 6 ' 120 ' 1680 ' 

= a n (l - x 2 ) + a,\ x — x 3 x 5 x 1 -■ 

\ 6 120 1680 



If we define 

,.,t . .. 1 3 1 5 I7 

y,(x) = l — x and y (x) = x — x x x 

Jl yi 6 120 1680 

then the general solution can be rewritten as y = a^y^x) + a l y 2 (x). 

Solved Problem 12.4 Determine whetherx = 0is a regular singular point 
of the differential equation 



Sx 2 y" + l0xy' + (x-l)y = 
Dividing by 8* 2 , we have 



J 1_ 

4x """ 8x Sx 2 



P(x) = — and Q(x) = — - n : 



Neither of these functions is defined at x = 0, so this point is a singular 
point. Furthermore, both 

xP(x) = ^ and x 2 Q(x) = ±(x-\) 
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are analytic every where: the first is a constant and the second a polyno- 
mial. Henee, both are analytic at x = 0, and this point is a regular singu- 
lar point. 

Solved Problem 12.5 Use the method of Frobenius to find one solution 
near x = of x 2 y" - xy' + y = . 

Here P(x) = -\/x and Q(x) = l/x 2 , so x = is a regular singular point 
and the method of Frobenius is appli cable. Substituting Equations 12.9 
through 12.11 into the left side of the differential equation, as given and 
combining coefficients of like powers of x, we obtain 

x x (X-\) 2 a +x x+l [X 2 a,] + --- 



+x x+n [(X + nf - 2{X + n) + \]a n +■■■ = () 



Thus, 



and, in general, 



(X-\) 2 a o =0 (12.19) 



[(X + n) 2 -2(X + n) + \]a n =0 (12.20) 



From Equation 12.19, the indicial equation is (X - 1 ) 2 = 0, which has roots 
X x = X 2 = 1. Substituting X = 1 into Equation 12.20 we obtain n 2 a n = 0, 
which implies that a n = 0, n > 1. Thus, y^x) = a^K. 

Solved Problem 12.6 Find the general solution near x = of ?>x 2 y" - 
xy' + y = . 

Here P(x) = -l/(3x) and Q(x) = l/(3x 2 ); henee, x = 0, is a regular sin- 
gular point and the method of Frobenius is applicable. Substituting Equa- 
tions 12.9 through 12.11 into the differential equation and simplifying, 
we have 

x Á [3X 2 -4X + l]a +x Á+] [3X 2 +2X]a 1 +- 
+x x+ "[3(X + n) 2 -4(X + n) + l]a n + --- = 

Dividing by x x and equating all coefficients to zero, we find 
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(3X 2 -U + \)a =0 (12.21) 



and 



[3(X + n) 2 -4(X + n) + l]a n =0 n>\ (12.22) 

From 12.21, we conclude that the indicial equation is 3X 2 - 4X + 1 = 0, 
which has roots X { = 1 and X = 1/3. Since X x — A 2 = 2/3, the solution is 
gi ven by Equations 12.12 and 12.13. Note that for either valué of X, Equa- 
tion 12.22 is satisfied by simply choosing a n = 0, n > 1. Thus, 

y 1 (x) = x 1 ^a n x n = a x 

n=0 

y 2 {x) = x m f j a n x n =a x m 

n=0 

and the general solution is 

y = c, v, (x) + c 2 y 2 (x) = k x x + k 2 x ] ' 3 

where k l =c l a Q and k 2 = c 2 a Q . 

Solved Problem 12.7 Use the method of Frobenius to find one solution 
near x = of x 2 y" + (x 2 - 2x)y' + 2y = 0. 

2 2 . . 

Here P(x) = 1 and Q{x) = —¿- ; so, x = is a regular singular 

xx 
point and the method of Frobenius is applicable. Substituting Equations 
12.9 through 12.11 into the left side of the differential equation, as giv- 
en, and combining coefficients of like powers of x, we obtain 

x x [(X 2 -3X + 2)a ] + x Á+] [(X 2 - X)a x + Xa ] + ■ ■ ■ 
+x x+ "{[(X + n) 2 -3(X + n) + 2]a n +(X + n-l)a n _ i } + --- = 

Dividing by x Á , factoring the coefficient of a n , and equating the coeffi- 
cient of each power of x to zero, we obtain 

(X 2 -3X + 2)a =0 (12.23) 
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and, in general, [(X + n) - 2] [(X + n)— X\a n + (X + n- V)a n _ x = 0, or, 

a n=- 7 J "n-x ("*1) (12.24) 

/l + n-2 

From 12.23, the indicial equation is X 2 - 3X + 2 = 0, which has roots X x 
= 2 and X 2 = 1 . Since Ij - A 2 = 1 , a positive integer, the solution is given 
by Equations 12.12 and 12.16. Substituting X = 2 into 12.24, we have a n 
= -( 1 /«)a . , from which we obtain 

a x = -a 

1 1 

a = — a, = — a 

2 2! 

1 11 1 

a-, = — a 7 = a n = a a 

3 3 2 32 , 3 , 

(-1/ 

and, in general, a t = a n . Thus, 

k k\ ° 

y x (x) = a x 2 £ — —x" = a x 2 e~ x 
n ni 

n=0 



Chapter 13 

Gamma 
and Bessel 
functions 



In this Chapter: 

• Gamma Function 

• Bessel Functions 

• Algébrate Operations on Infinite 
Series 

• Solved Problems 

Gamma Function 

The gamma function, T(p), is defined for any positive real number/? by 

T(p) = jx p - X e- X dx (13.1) 

o 
Consequently, T(l) = 1 and for any positive real number/?, 

r(p + \) = pT(p) (13.2) 

Furthermore, when p = n, a positive integer, 

T(n + l) = n\ (13.3) 

98 
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Thus, the gamma function (which is defined on all positive real numbers) 
is an extensión of the factorial function (which is defined only on the non- 
negative integers). 

Equation 13.2 may be rewritten as 

T(p) = -r(p + \) (13.4) 

P 

which defines the gamma function iteratively for all nonintegral negative 
valúes of p. T(0) remains undefined, because 

hm i(p)= lim — = oo 

p-»0 + p->0 + p 

and 

1- TV ^ I' r 0> + D 

hm 1 (/?) = hm = —°° 

p->0~ p-»o~ p 

It then folio ws from Equation 13.4 that T(p) is undefined for negative in- 
teger valúes of p. 

Table 13.1 lists valúes of the gamma function in the interval 1 <p < 
2. These tabular valúes are used with Equations 13.2 and 13.4 to genér- 
ate valúes of T(p) in other intervals. 

Bessel Functions 

Let p represent any real number. The Bessel function of the first kind of 
orcler p, J (x), is 

JJx) = T „} ' (13.5) 

The function J pe) is a solution near the regular singular point x = of 
Bessel's differential equation oforderp: 

x 2 y" + xy' + (x 2 -p 2 )y = (13.6) 

In fact, / (x) is that solution of Equation 13.6 guaranteed by Theorem 
12.2. 
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X 


r(A-) 


X 


r(x) 


X 


r(.i-) 


1.00 


1.0000 0000 


1.34 


0.8922 1551 


1.67 


0.9032 9650 


1.01 


0.9943 2585 


1.35 


0.89115144 


1.68 


0.9050 0103 


1.02 


0.9888 4420 


1.36 


0.8901 8453 


1.69 


0.9067 8182 


1.03 


0.9835 4995 


1.37 


0.8893 1351 


1.70 


0.9086 3873 


1.04 


0.9784 3820 


1.38 


0.8885 3715 


1.71 


0.9105 7168 


1.05 


0.9735 0427 


1.39 


0.8878 5429 


1.72 


0.9125 8058 


1.06 


0.9687 4365 


1.40 


0.8872 6382 


1.73 


0.9146 6537 


1.07 


0.9641 5204 


1.41 


0.8867 6466 


1.74 


0.9168 2603 


1.08 


0.9597 2531 


1.42 


0.8863 5579 


1.75 


0.9190 6253 


1.09 


0.9554 5949 


1.43 


0.8860 3624 


1.76 


0.9213 7488 


1.10 


0.9513 5077 


1.44 


0.8858 0506 


1.77 


0.9237 6313 


1.11 


0.9473 9550 


1.45 


0.8856 6138 


1.78 


0.9262 2731 


1.12 


0.9435 9019 


1.46 


0.8856 0434 


1.79 


0.9287 6749 


1.13 


0.9399 3145 


1.47 


0.8856 3312 


1.80 


0.9313 8377 


1.14 


0.9364 1607 


1.48 


0.8857 4696 


1.81 


0.9340 7626 


1.15 


0.9330 4093 


1.49 


0.8859 4513 


1.82 


0.9368 4508 


1.16 


0.9298 0307 


1.50 


0.8862 2693 


1.83 


0.9396 9040 


1.17 


0.9266 9961 


1.51 


0.8865 9169 


1.84 


0.9426 1236 


1.18 


0.9237 2781 


1.52 


0.8870 3878 


1.85 


0.9456 1118 


1.19 


0.9208 8504 


1.53 


0.8875 6763 


1.86 


0.9486 8704 


1.20 


0.9181 6874 


1.54 


0.8881 7766 


1.87 


0.9518 4019 


1.21 


0.9155 7649 


1.55 


0.8888 6835 


1.88 


0.9550 7085 


1.22 


0.9131 0595 


1.56 


0.8896 3920 


1.89 


0.9583 7931 


1.23 


0.9107 5486 


1.57 


0.8904 8975 


1.90 


0.9617 6583 


1.24 


0.9085 2106 


1.58 


0.8914 1955 


1.91 


0.9652 3073 


1.25 


0.9064 0248 


1.59 


0.8924 2821 


1.92 


0.9687 7431 


1.26 


0.9043 9712 


1.60 


0.8935 1535 


1.93 


0.9723 9692 


1.27 


0.9025 0306 


1.61 


0.8946 8061 


1.94 


0.9760 9891 


1.28 


0.9007 1848 


1.62 


0.8959 2367 


1.95 


0.9798 8065 


1.29 


0.8990 4159 


1.63 


0.8972 4423 


1.96 


0.9837 4254 


1.30 


0.8974 7070 


1.64 


0.8986 4203 


1.97 


0.9876 8498 


1.31 


0.8960 0418 


1.65 


0.9001 1682 


1.98 


0.9917 0841 


1.32 


0.8946 4046 


1.66 


0.9016 6837 


1.99 


0.9958 1326 


1.33 


0.8933 7805 











Table 13.1 
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Algebraic Operations on Infinite Series 

Changing the dummy índex. The dummy index in an infinite series can 
be changed at will without altering the series. For example, 



1 777777 = 1777777= I 



1 11111 

- = — + — + — + — + — + • 



tj(* + l)! fo( n + i y- ¿ToO> + 1 ) ! 1! 2! 3! 4! 5! 

V l 
Change of variables. Consider the infinite series /_, ~ — ~- If we make 

. £(* + 1 ) 1 

the change of variables y = k + 1, or k =j - 1 , then 

É— =1- 

Note that a change of variables generally changes the limits on the sum- 
mation. For instance, ify = k + 1, it follows that y = 1 when k = 0,y = <*>, 
when k = °°, and, as k runs from to °o,y* runs from 1 to <*>. 

The two operations given above are often used in concert. For ex- 
ample, 

f_^_ = y_}_ = y_}_ 

Here, the second series results from the change of variables j = k + 2 in 
the first series, while the third series is the result of simply changing the 
dummy index in the second series from y to k. Note that all three series 
equal 

1111 

— + — + — + — + ••• = £-1 
1! 2! 3! 4! 



Solved Problems 

Solved Problem 13.1 Pro ve that T{p + 1) =pT(p), p>0. 
Using Equation 13.1 and integration by parts, we have 
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T(p + \)= ¡x (p+1) - l e- x dx= lim [x p e- x dx 



= lim 



-x p e- x \ + 



jpx 



P- X o~ x 



e~ x dx 



= üm(-r p e- r +0) + p\x p - 1 e- x dx = pT(p) 



The result lim r p e ' " = is easily obtained by first writing r p e '" as r p le r 
and using L'Hópital's rule. 

Solved Problem 13.2 Use the method of Frobenius to find one solution 
of Bessel's equation of order p: 

x 2 y" + xy' + (x 2 -p 2 )y = 

Substituting Equations 12.9 through 12. 1 1 into Bessel's equation and 
simplirying, we find that 

x\X 2 -p 2 )a +x x+x [{X + \) 2 -p 2 ]a x + 
+x Á+2 {[(X + 2) 2 -p 2 ]a 2 +a } + - 

+x Á+ "{[(X + n) 2 -p 2 ]a n +a n _ 2 } + - = 



Thus, 



(X 2 -p 2 )a o =0 U + \f-p 2 ]a x =0 (13.7) 



and, in general, [(X + n) 2 - p 2 ]a n + a n _ 2 = 0, or 

1 



a„=- 



(X + nf-p 2 



a n _ 2 (n>2) 



(13.8) 



The indicial equation is X 2 - p 2 = 0, which has roots X x = p and X 2 = -p 
(p nonnegative). Substituting X = p into 13.7 and 13.8 and simplirying, 



we find that a , = and 



a„ = — 



J 

n(2p + n) 



a n -2 (n>2) 
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Henee, = a, =a 3 =a 5 =a 7 =■■■ and 



a n = 



-1 



2 2 2 l!(/7 + l) ° 



a. =-- 



-a, =• 



4 2 2 2(p + 2) 2 2 4 2\(p + 2)(p + l) ° 



-1 



a* =-- 



-a A =■ 



6 2 2 3(/? + 3) 4 2 6 3!(/7 + 3)(/7 + 2)(/p + 1) ° 



and, in general, 



a,, =■ 



(-lf 



2k 2 2k k\(p + k)(p + k-l)---(p + 2)(p + \)" í 



a (k>\) 



y l (x) = x x J j a n x n = x l 



•X 



«0+2A* 
k = i 



2k 



Thus, 



= a a x' 



l+Z: 



(-1)** 2 * 



, =1 '2 2 ^!(/7 + A:)(/7 + A:-l)-(/7 + 2)(p + l) 

It is customary to choose the arbitrary constant a Q as a = 

Then bringing a^^ inside the bracket 
finally using Problem 13.1, we obtain 



1 



" u 2 p r(/?+i) 

Then bringing a^^ inside the brackets and summation, combining, and 



2T(/? + l) ~2 2A+/, A:!r(/? + A: + l) 






*= 



w 2 2/r +Pit!r(/7+it + 1) 
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Direction Fields 

Graphical methods produce plots of solutions to first-order differential 
equations of the form 

y' = f(x,y) (14.1) 

where the derivative appears only on the left side of the equation. 

Example 14.1 (a) For the problem y' = — y + x + 2, we have/(x, y) = —y 
+ x + 2. (b) For the problem y' =y+ 1, we nave /(a;, y) = y 2 + 1 . (c) For 
the problem y' = 3, we have/(jt, y) = 3. Observe that in a particular prob- 
lem, f(x, y) may be independent of x, of y, or of x and y. 

Equation 14.1 defines the slope of the solution curve y(x) at any point 
(x, y) in the plañe. A line element is a short line segment that begins at the 
point (x, y) and has a slope specified by 14.1; it represents an approxi- 
mation to the solution curve through that point. A collection of line ele- 
ments is a direction fie Id. Note that direction fields are tedious to draw by 
hand, but they can be easily created by many computer algebra systems. 
The graphs of solutions to 14. 1 are generated from direction fields by 
drawing curves that pass through the points at which line elements are 
drawn and also are tangent to those line elements. 



ik 



Note! 



A direction field can be used to graphically deter- 
mine the behavior of the solution. 



If the left side of Equation 14. 1 is set equal to a constant, the graph of the 
resulting equation is called an isocline. Different constants define differ- 
ent isoclines, and each has the property that all elements emanating from 
points on that isocline have the same slope, a slope equal to the constant 
that generated the isocline. When they are simple to draw, isoclines yield 
many line elements at once which are useful for constructing direction 
fields. 
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Euler's Method 



If an initial condition of the form 



y( x o)=y 



(14.2) 



is also specified, then the only solution curve of Equation 14.1 of inter- 
est is the one that passes through the initial point (x Q , y ). 

To obtain a graphical approximation to 
the solution curve of Equations 14.1 and 
14.2, begin by constructing a line element at 
the initial point (x Q , v ) and then continuing 
it for a short distance. Denote the terminal 
point of this line element as (x l ,y l ). Then 
construct a second line element at(x v y l ) and 
continué it a short distance. Denote the ter- 
minal point of this second line element as (x 2 , 
y 2 ). Follow with a third line element con- 
structed at (x 2 , y 2 ) and continué it a short distance. The process proceeds 
iteratively and concludes when enough of the solution curve has been 
drawn to meet the needs of those concerned with the problem. 

If the difference between successi ve x valúes are equal, that is, if for 
a specified constant h, h=x l -x Q = x 2 ~x l =x 3 -x 2 = ■■■, then the graph- 
ical method given above for a first-order initial- valué problem is known 
as Euler's method. It satisfies the formula 




yn + i=yn +h f(x n ,y n ) 

for n = 1,2,3,.... This formula is often written as 

yn +i =yn +h y'n 

where 

as required by Equation 14.1. 



(14.3) 



(14.4) 



(14.5) 
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Stability 

The constant h in Equations 14.3 and 14.4 is called the step-size, and its 
valué is arbitrary. In general, the smaller the step-size, the more accurate 
the approximate solution becomes at the price of more work to obtain that 
solution. Thus, the final choice of h may be a compromise between ac- 
curacy and effort. If h is chosen too large, then the approximate solution 
may not resemble the real solution at all, a condition known as numeri- 
cal instability. To avoid numerical instability, Euler's method is repeat- 
ed, each time with a step-size one-half its previous valué, until two suc- 
cessive approximations are cióse enough to each other to satisfy the needs 
of the solver. 

General Remarks Regarding 
Numerical Methods 

A numerical method for solving an initial- valué problem is a procedure 
that produces approximate solutions at particular points using only the 
operations of addition, subtraction, multiplication, división, and func- 
tional evaluations. In this chapter, we consider only first-order initial- val- 
ué problems of the form 

Y = f{x,y)\ y{x ) = y (14.6) 

Generalizations to higher-order problems are given later in this chapter. 
Each numerical method will produce approximate solutions at the points 
x Q ^c l ^c 2 ,..., where the difference between any two successi ve x- valúes is 
a constant step-size h: that is x , . - x = h(n = 0,1,2,...). Remarks made 

f 5 «+I ti V 5 5 5/ 

previously in this chapter on the step-size remain valid for all the nu- 
merical methods presented. 

The approximate solution at x n will be designated by y(x n ), or sim- 
ply y . The true solution at x will be denoted by either Y(x ) or Y . Note 

r J J n n J v n' n 

that once y n is known, Equation 14.6 can be used to obtain y' n as 

y'„ = f(x n >y n ) (14.7) 

The simplest numerical method is Euler's method, described earlier 
in this chapter. 
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A predictor-corrector method is a set of two equations for y n¥V The 
first equation, called the predictor, is used to predict (obtain a first ap- 
proximation to) y n+ {, the second equation, called the corrector, is then 
used to obtain a corrected valué (second approximation) to y n+v In gen- 
eral, the corrector depends on the predicted valué. 

Modified Euler's Method 

This is a simple predictor-corrector method that uses Euler's method as 
the predictor and then uses the average valué of / at both the left and right 
end points of the interval \x n , x rH _ 1 \ (n = 0,1,2,...) as the slope of the line 
element approximation to the solution over that interval. The resulting 
equations are: 

predictor: y n+l =y n +hy' n 
corrector: y n+1 =y n + -{y' n+x + y' n ) 

For notational convenience, we designate the predicted valué of y n+l by 
py n+l - It then follows from Equation 14.7 that 

py'n + i=f(Xn + i>py n+ i) (14-8) 

The modified Euler's method becomes 

predictor: py n+l =y n +hy' n 

h, , ,, (14.9) 

corrector: y n+x = y n + - (py n+] + y n ) 

Runge-Kutta Method 

>Wi =y n +^(k l +2k 2 +2k,+k 4 ) (14.10) 



where 



K=hf{x n ,y n ) 

k 2 =hf(x n +-h,y n +-k l ) 

k3=hf(x n +-h,y n +-k 2 ) 
k 4 = hf(x n + h,y n +k 3 ) 
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This is not a predictor-corrector method. 

Adams-Bashforth-Moulton Method 

predictor: py n+] = y n + ^-(55/ n -#/_, +37)/_ 2 -9^_ 3 ) 

/T (14.11) 

corrector: y n+l = y n +—(9py' n+1 +I9y' n - 5y' n _ x + y' n _ 2 ) 

Milne's Method 

4/i 
predictor: /?y„ +I = y n _ 3 + — {2y' n - y' n _ x + 2y' n _ 2 ) 

corrector: y n+l = y n _ x +-(py n+l + 4y n + >■„_, ) 

Starting Valúes 

The Adams-Bashforth-Moulton method and Milne's method require in- 
formaron at y ,y v y 2 > an d ^ t0 start. The first of these valúes is given by 
the initial condition in Equation 14.6. The other three starting valúes are 
gotten by the Runge-Kutta method. 

Order of a Numerical Method 

A numerical method is of order n, where n is a positi ve 
integer, if the method is exact for polynomials of de- 
gree n or less. In other words, if the true solution of an 
initial- valué problem is a polynomial of degree n or 
less, then the approximate solution and the true solu- 
tion will be identical for a method of order n. 

In general, the higher the order, the more accurate the method. 
Euler's method, Equation 14.4, is of order one, the modified Euler's 
method, Equation 14.9, is of order two, while the other three, Equations 
14.10 through 14.12, are fourth-order methods. 
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Numerical Methods for Systems 

First-Order Systems 

Numerical methods for solving first-order initial- valué problems, includ- 
ing all of those described previously in this chapter, are easily extended 
to a system of first-order initial- valué problems. These methods are also 
applicable to most higher-order initial-value problems, in particular 
those that can be transformed to a system of first-order differential equa- 
tions by the reduction process described in Chapter Eleven. 
Standard form for a system of two equations is 

y' = f(x,y,z) 

z' = g(x,y,z); (14.13) 

y(x ) = y , z(x ) = z 

We note that, with y' = f(x,y,z) = Z, with system 14.13 represents the 
second-order initial-value problem 

y" = g(x,y,yj, y(x ) = y , y'(x ) = z 

Standard form for a system of three equations is 

Y = f(x, y,z,w) 
z' = g(x,y,z,w) 
w' = r{x,y,z,w)\ (14.14) 

y(x ) = yo> z(x ) = z , w(x )=w 

If, in such a system, f(x, y, z, w) = z, and g(x, y, z, w) = w, then system 
14. 14 represents the third-order initial-value problem 

y"' = r(x,y,z,w); y(x ) = y , y'(x ) = Zo, y"(x ) = w 

The formulas that follow are for systems of two equations in stan- 
dard form (14.13). Generalizations to systems of three equations in stan- 
dard form (14.14) or systems with four or more equations are straight- 
forward. 
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Euler's Method 



y n+ \=yn+ h yn 

Z n +\ = Z n + "^« 



Runge-Kutta Method 



where 



(14.15) 



y n+ i=y n +T( k i +2k 2+^+k 4 ) 

1 (14.16) 

z n+x =z n +-d x +2l 2 +2l 3 +l 4 ) 
6 



k x=hf(x n ,y n ,z n ) 

h= h 8( X n>y n >Zn) 

k 2 =hf{x n + \h,y n + \k li z n + \h) 
k=hg(x n + ih,y n + ^k l ,z n +il 1 ) 
^ = ¥(x n +\h,y n + ^k 2 ,z n + \l 2 ) 
h=hg(x n + ±h iyn + \k 2 ,z n + il 2 ) 
k 4 =hf(x n + h,y n +k 3 ,z n + l 3 ) 
U=hg(x n +h,y n +k 3 ,z l7 +l 3 ) 



Adams-Bashforth-Moulton Method 



h 
predictors: py n+l =y n + — (55y' n - 59y' n _ l + 37^_ 2 - 9y' n _ 3 ) 

h 

PZ n+ x = z n +^( 55 Zn ~ 59z ; ' 7 _, + 374_ 2 - 9^3 ) 

correctors: y n+1 = y n +—(9py' n+l + \9y' n -5y' n _ ] +y' n _ 2 ) ( l4 ' 7) 

h 
z n+ i =z n +—(9pz' n+1 +19z' n -5z' n .. i +z' n - 2 ) 

Corresponding derivatives are calculated from system 14.13. In particu- 
lar, 
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^w+1 = 8\ X n+\*yn+\*Zn+l) 



(14.18) 



The derivatives associated with the predicted valúes are obtained simi- 
larly, by replacing y and z in Equation 14.18 with py and pz, respective- 
ly. As in the previous section, four sets of starting valúes are required for 
the Adams-Bashforth-Moulton method. The first set comes directly from 
the initial conditions; the other three sets are obtained from the Runge- 
Kutta method. 

Solved Problems 

Solved Problem 14.1 Draw two solution curves to the differential equa- 
tion y = 2y — x. 

A direction field for this equation is given by Figure 14-1. Two solution 
curves are also shown, one that passes through the point (0,0) and a sec- 
ond that passes through the point (0,2). Observe that each solution curve 
folio ws the flow of the line elements in the direction field. 
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Figure 14-1 
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Solved Problem 14.2 Find y(l) for y r = y — x; y(0) = 2, using Euler's 
method with h = \. 

For this problem, x Q = 0, y = 2, and/(jc, y) = y - x; so Equation 14.5 be- 
comes y' n = y n — x n . Because h = \, 

x x = x + h = j x 2 = x l +h = ±- x 3 =x 2 + h = j x 4 = x 3 + h = l 

Using Equation 14.4 with n = 0,1,2,3 successively, we now compute the 
corresponding y- valúes. 

n = 0: yi=y +fy 

But^=/(x , y ) = y -x = 2 - = 2 
Henee, tt=2 + l(2) = f 

n=1: ) 7 2 = ) ; i +/í /l 

Buty i =f(x v y 1 ) = y 1 -x 1 = | -■}■ = ■£ 
Henee, y 2 = 1 + L ( 1 ) = £ 

n = 2: y 3 = y 2 + hy 2 

_ 49 1 _ 41 



Buty 2 =f(x 2 ,y 2 ) = y 2 -x 2 = %-± = % 

V =49+1(11)= 237 
"3 16 4M6' 64 



Henee v = ü + ±fll) = 237 



n = 3: y 4 = y 3 + h/ 3 

Rut V = f(x V1=V -X =237_2-189 
Upnrp v - 237 , 1 / 1 89 \ _ 1137 

Henee, y 4 - _ + -(_)_^- 



Thus, 



^l)*^^* 4.441 



Note that the true solution is Y(x) = ex + x + 1 , so that F(l) « 4.7 1 8. 
If we plot (x , v ) for « = 0,1,2,3, and 4, and then connect successive 
points with straight line segments, as done in Figure 14-2, we have an ap- 
proximation to the solution curve on [0,1] for this initial- valué problem. 
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*-x 



Figure 14-2 



Chapter 15 

B ound ar y- Valué 
Problems and 
Fourier Series 



In This Chapter: 

• Second-Order Boundary- Valué 
Problems 

• Eigenvalue Problems 

• Sturm-Liouville Problems 

• Eigenfunction Expansions 

• Fourier Sine Series 

• Fourier Cosine Series 

• Solved Problems 

Second-Order Boundary-Value Problems 

Standard Form 

A boundary- valué problem in standard form consists of the second-order 
linear differential equation 

y" + P{x)y f + Q{x)y = <j>{x) (15.1) 
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and the boundary conditions 

a l y(a) + /3 l y'(a) = y l 
a 2 y(b) + p 2 y'(b) = y 2 

where P(x), Q(x), and (¡)(x) are continuous in [a, b] and a v a 2 , P V P 2 , y v 
and y 2 are all real constants. Furthermore, it is assumed that a l and/?, are 
not both zero, and also that a and fi 2 are not both zero. 

The boundary- valué problem is said to be homogeneous if both the 
differential equation and the boundary conditions are homogeneous (i.e., 
<p(x) = and y x = y 2 = 0). Otherwise the problem is nonhomogeneous. 
Thus a homogeneous boundary- valué problem has the form 

y" + P(x)y' + Q(x)y = 0; 

aí y(a) + p iy '(a) = (15.3) 

a 2 y(b) + p 2 y'(b) = 

A somewhat more general homogeneous boundary- valué problem than 
15.3 is one where the coefficients P(x) and Q(x) also depend on an arbi- 
trary constant X. Such a problem has the form 

y" + P(x,l)y' + Q{x,X)y = 0; 

a l y(a)+p i y(a) = (15.4) 

a 2 y(b)+p 2 y'(b) = 

Both 15.3 and 15.4 always admit the trivial solution y(x) = 0. 

Solutions 

A boundary- valué problem is solved by first obtaining the general solu- 
tion to the differential equation, using any of the appropriate methods pre- 
sented heretofore, and then applying the boundary conditions to evalúate 
the arbitrary constants. 

Theorem 15. 1. Let y x (x) and y 2 (x) be two linearly independent solutions 
of 

y" + P(x)y' + Q(x)y = 
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Nontrivial solutions (i.e., solutions not identically equal to zero) to the 
homogeneous boundary- valué problem 15.3 exist if and only if the de- 
terminant 



a l y l (a) + P x y[(a) a l y 2 (a) + Pj 2 (a) 
a 2 y l (b) + P 2 y[(b) a 2 y 2 (b) + p 2 y' 2 (b) 



(15.5) 



equals zero. 



Theorem 15.2. The nonhomogeneous boundary- valué problem defined 
by 15.1 and 15.2 has a unique solution if and only if the associated ho- 
mogeneous problem 15.3 has only the trivial solution. 

In other words, a nonhomogeneous problem has a unique solution 
when and only when the associated homogeneous problem has a unique 
solution. 

Eigenvalue Problems 

When applied to the boundary- valué problem 15.4, Theorem 15.1 shows 
that nontrivial solutions may exist for certain valúes of X but not for oth- 
er valúes of X. 



You Need to Know 



/ 



Those valúes of X for which nontrivial solutions do 
exist are called eigenvalues; the corresponding 
nontrivial solutions are called eigenfunctions. 



Sturm-üouville Problems 

A second-order Sturm-Liouville problem is a homogeneous boundary- 
value problem of the form 

\j>(x)/Y + q(*)y + Xw(x)y = 0; (15.6) 
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a x y{a) + p x y'(a) = 

a 2 y(b) + p 2 y'(b) = 

where p(x), p'(x), q(x), and w(x) are continuous on [a, b], and both p{x) 
and w(x) are positive on [a, b\. 

Equation 15.6 can be written in standard form 15.4 by dividing 
through by p(x). Form 15.6, when attainable, is preferred, because Sturm- 
Liouville problems have desirable features not shared by more general 
eigenvalue problems. The second-order differential equation 

a 2 (x)y" + a x {x)y' + a (x)y + Xr{x)y = (15.8) 

where a 2 (x) does not vanish on [a, b], is equivalent to Equation 15.6 if 
and only if a' 2 (x) = a { (x). This condition can always be forced by multi- 
plying Equation 15.8 by a suitable factor. 

Properties of Sturm-Liouville Problems 

Property 15.1. The eigenvalues of a Sturm-Liouville problem are real 
and nonnegative. 

Property 15.2. The eigenvalues of a Sturm-Liouville problem can be 
arranged to form a strictly increasing infinite sequence; that is, < X x < 
X- < X- < • • • . Furthermore, X — > °° as n — > ~. 

2 3 ' n 

Property 15.3. For each eigenvalue of a Sturm-Liouville problem, there 
exists one and only one linearly independent eigenfunction. 

[By Property 15.3 there corresponds to each eigenvalue X n a unique 
eigenfunction where the coefficient on the lead term is one; we denote 
this eigenfunction by e n (x).] 

Property 15.4. The set of eigenfunctions {e^x), e 2 (x),...} of a Sturm- 
Liouville problem satisfies the relation 



\w(x)e n (x)e m (x)dx = (15.9) 

a 

for n^m, where w(x) is given in Equation 15.6. 
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Eigenfunction Expansions 

A wide class of functions can be represented by infinite series of eigen- 
functions of a Sturm-Liouville problem. 

Definition: A function f{x) is piecewise continuous on the open interval 
a <x<b if (l)/(x) is continuous everywhere in a < x < b with the possi- 
ble exception of at most aftnite number of points x v x 2 ,..., x n , and (2) at 
these points of discontinuity, the right- and left-hand limits of f(x) re- 
spectively lim f(x) and lim f(x), existí; = 1,2,..., n). 

X— >Xj X^tXj 

x>x x<x 

(Note that a continuous function is piecewise continuous.) 

Definition: A function /(x) is piecewise continuous on the closed inter- 
val a<x<bif(l)itis piecewise continuous on the open interval a<x<b, 
(2) the right-hand limit of f(x) exists at x = a and (3) the left-hand limit 
of f(x) exists at x = b. 

Definition: A function /(x) is piecewise smooth on [a, b] if both/(x) and 
f(x) are piecewise continuous on [a,b]. 

Theorem 15.3. lff(x) is piecewise smooth on [a, b] and if [e n {x)} is the 
set of all eigenfunctions of a Sturm-Liouville problem (see Property 
15.3), then 



where 



/W=Iv»W (15.10) 



jw(x)f(x)e n (x)dx 
c n =— b (15.11) 



\w(x)el{x)dx 



The representation 15.10 is valid at all points in the open interval {a, b) 
where /(x) is continuous. The function w(x) in 15.11 is given by Equa- 
tion 15.6. 
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Because different Sturm-Liouville problems usually genérate differ- 
ent sets of eigenfunctions, a given piecewise smooth function will have 
many expansions of the form 15.10. The basic features of all such ex- 
pansions are exhibited by the trigonometric series discussed below. 

Fourier Sine Series 

The eigenfunctions of the Sturm-Liouville problem y" + Xy = 0; v(0) = 0, 
y(L) = 0, where L is a real positive number, are ejx) = ún{nnxlL) (n = 
1,2,3,...). Substituting these functions into 15.10, we obtain 

/(*) = ¿ C;? sin^ (15.12) 

For this Sturm-Liouville problem, w(x) =1, a = 0, and b = L; so that 

b L 

\ w(x)e^(x)dx = \ sin 2 dx = — 

a 

and 15.11 becomes 

? L 
c n =-\f{x)ún — dx (15.13) 

L o L 

The expansión 15. 12 with coefficients given by 15. 13 is the Fourier sine 
series for/(x) on (0, L). 

Fourier Cosine Series 

The eigenfunctions of the Sturm-Liouville problem y" + Xy = 0; /(0) = 
0, v(L) = 0, where L is a real positive number, are e Q (x) = 1 and e n (x) = 
cos(rmx/L) (n = 1,2,3,...). Here X = is an eigenvalue with correspond- 
ing eigenfunction e Q (x) = 1. Substituting these functions into 15.10, 
where because of the additional eigenfunction e Q (x) the summation now 
begins at n = 0, we obtain 

f(x) = c +^c„cos?j- (15.14) 

n=\ 
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For this Sturm-Liouville problem, w(x) = 1, a = 0, and b = L; so that 



f w(x)eQ (x)dx = f dx = L 



and 



f w(x)e„ (x)dx = f eos 2 dx = — 



Thus 15.11 becomes 

L ^ L 



c =j]f(x)dx c n =jjf(x)cos^dx (15.15) 

(n = 1,2,...) 



o o 



The expansión 15.14 with coefficients given by 15.15 is the Fourier co- 
sine series foxf(x) on (0, L). 

Solved Problems 

Solved Problem 15.1 Solve /' + 2/ - 3;y = 9x; y(0) = 1, /(l) = 2. 

This is a nonhomogeneous boundary- valué problem of forms 15.1 
and 15.2, where (p(x) = 9x,y l = \, and y = 2. The general solution to the 
homogeneous differential equation is v = c x e~ 3x + c 2 e*. If we apply ho- 
mogeneous boundary conditions, we find that c x = c 2 = 0; henee the so- 
lution to the homogeneous problem is y = 0. Since the associated homo- 
geneous problem has only the trivial solution, it follows from Theorem 
1 5.2 that the given problem has a unique solution. Solving the differential 
equation by the method of Chapter Six, we obtain 

y = c { e~ 3x +c 2 e x -3x-2 

Applying the boundary conditions, we find 
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c l +c 2 -2 = 1 -3c x e +c 2 e-3 = 2 



whence 



c, = 



3^-5 



Co =■ 



5 + 9e~ 



1 - ->-3 2 — „ -3 

¿ + 3é? e + 3<? 



Finally, 



3 , = - 



(3e-5y _3v +(5 + 9^ 3 )^ 



e + 3e 



-3x-2 



Solved Problem 15.2 Find the eigenvalues and eigenfunctions of 

/'-4¿/ + 4A 2 ;y = 0; y(0) = 0, >'(1) + /(1) = 

The coefficients of the given differential equation are constants (with 
respect to x); henee, the general solution can be found by use of the char- 
acteristic equation. We write the characteristic equation in terms of the 
variable m, since X now has another meaning. Thus we have m 2 - 4Xm + 
4X 2 = 0, which has the double root m = 2X; the solution to the differential 
equation is y = c x e 2Xx + c^ce 2Xx . Applying the boundary conditions and 
simplifying, we obtain 

c 1 = c l (l + 2X) + c 2 (2 + 2X) = 



It now folio ws that c l = and either c 2 = or X = — 1. The choice 
c 2 = results in the trivial solution y = 0; the choice X = -1 results in the 
nontrivial solution y = c^e 2 *, c 2 arbitrary. Thus, the boundary- valué 
problem has the eigen valué X = -1 and the eigenfunction y = c^xe 2 *. 

Í0 x<2 
Solved Problem 15.3 Find a Fouriersine series for f(x) = < on 

(0,3). [2 x>2 



Using Equation 15.13 with L = 3, we obtain 



z f ti \ ■ mx i 
c„=-J/O0sm — " v 

3 o 



3 
nnx 



= — í(0)sin : -^-í/v + 



4J (2)si 



nnx , 
sin dx 



= + ¿ 
3 



3 nnx 
eos 

nn 3 



2 

3J V 3 

2 
*=3 4 



2«7T 

eos coswr 

3 
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Thus Equation 15.12 becomes 



2nn „ 

eos (-1) 

3 



sin- 



nitx 



„ . 2n 1 4n 1 6tü . 

Hirthermore, eos — = — , eos — = — , eos — = 1,... 

3 2 3 2 3 



Henee, 



... 4(1 . Toe 3 . 2xx 2 . ?>kx 

f(x) = —\ —sin sin h— sin 

J k\2 3 4 3 3 3 



Since/(x) is piecewise smooth on [0,3] and continuous everywhere in 
(0,3) except at x = 2, it follows from Theorem 15.3 that this equality is 
valid everywhere in (0,3) except at x = 2. 



Appendix 

Laplace 
Transforms 





A*) 


F( S ) = X{/(*)} 


1. 


i 


- (J>0) 

5 


2. 


X 


4 (* >o ) 
5 


3. 


x nA (n = 1,2,...) 


( - 1)! C>0) 


4. 


V^ 


~4ñs~ i!1 (s>Q) 


5. 


i/Vjc 


4ñs~ xl1 (s>0) 


6. 


y M/2 (« = i,2,...) 


(l)(3)(5)--(2«-l)^„. B _ 1/2 

r 

(s>0) 


7. 


e <u 


(J>«) 


8. 


sin ax 


fl í"c ^ ni 


2 2 l J " U J 


9. 


eos ax 


5 (?^(Y\ 


2 2 V' V) 
S + fl 


10. 


sinh <7X 


^ / 1 IV 


2 2 («s - |ai; 
í -fl 
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/(*) 


F(s) = £{/(*)} 


ll. 


cosh ax 


* /" I |V 


2 2 (-S " W) 

s -a 


12 


x sin ax 


las 


(^ + « 2 ) 2 ( " 0) 


13. 


x eos ax 


2 2 

5 " a (--O) 


(,'+„*)' (J ' 0) 


14. 


*"V V (« = 1,2, ...) 


í- 1 » (.,>«) 
(*-«)" 


15. 


e sin ax 


*^ 1N 


(*-*)> +** °" 6) 


16. 


e' x eos ax 


í — ¿ , 


(,-6) 2 + « 2 (S ^ 


17. 


sin ax - a* eos ax 


2Cl3 (r^G\ 


(s 2 +* 2 ) 2 ( ) 


18. 


-e~ x/a 
a 


1 
1 + £75 


19. 




1 


s(j - ¿z) 


20. 


l- e - Wfl 


1 


s(\ + ¿w) 


21. 


<2 


1 


(\ + as) 2 
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A*) 


F(s) = X{f[x)} 


22. 




ax hx 

e -e 




1 


a-b 


(s - a)(s -b) 


23. 




e -xla _ e -xlb 




1 


a-b 


(l + as)(l + bs) 


24. 


{\+ax)e ax 


s 


(s-af 


25. 


a 


s 


(\ + as) 2 


26. 


ae ax -be bx 
a -b 


s 


(s - a)(s - b) 


27. 


i 


ie- xlb -be~ x!a 




s 


ab(a-b) 


(\ + as)(\ + bs) 


28. 


a 


1 


s 2 (s-a) 


29. 


sin 2 ax 


2a 2 


s(s 2 +4a 2 ) 


30. 


sinh 2 ax 


la 2 


s(s 2 -4a 2 ) 


31. 


1 

71 


. ax . ax 
cosh— p^sin— j= - 

V2 V2 

. , ax ax 
sinh— prcos— j= 

K V2 42) 




a 3 


s 4 +a 4 


32. 


. ax . . ax 
sin sinn 


a 2 s 




4% 4í 




s 4 +a 4 
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Ax) 


F(*)=X{/(x)} 


33. 


i 


ax . , ax 
eos— =• sinh— =■ + 

V2 V2 

. ax , ax 
sin — — cosh — ¡= 

v. V2 V2j 




a. 2 


5 4 +a 4 


34. 


ax . ax 

eos — -=r cosh —-= 

4i 4i 


s 3 


4 4 

j +a 


35. 


— (sinh ax - sin ax) 


3 

a 


s 4 -a 4 


36. 


— (cosh ax - eos ax) 


a 2 s 


s 4 -a 4 


37. 


— (sinh ax + sin ax) 


as 2 


s 4 -a 4 


38. 


— (cosh ax + eos a*) 


, 3 


4 4 

s — a 


39. 


sin ax sinh ax 


2a s 

s 4 +4a 4 


40. 


eos ax sinh ax 


a{s 2 -2a 2 ) 
s 4 +4a 4 


41. 


sin ax cosh ax 


a(s 2 +2a 2 ) 
s 4 +4a 4 


42. 


eos ax cosh ax 


s 3 


s 4 +4a 4 


43. 


1 ^ 


as 2 


2 






(s 2 +a 2 ) 2 
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A*) 


F(s) = £{f(x)} 


44. 


ax . 
eos ax sin ax 

2 


s 2 


(s 2 +a 2 ) 2 


45. 


— (ax cosh ax - sinh ax) 
2 


a 3 


U 2 -« 2 ) 2 


46. 


— sinh ax 

2 


as 


(* 2 -* 2 ) 2 


47. 


~ (sinh ax + ax cosh ax) 


as 


í 2 2\2 

-a ) 


48. 


cosh ax + — sinh ax 

2 


5 3 


(* 2 -« 2 ) 2 


49. 


a sin bx — b sin ax 


a¿> 


a 2 -6 2 


(í 2 +« 2 )(j 2 +6 2 ) 


50. 


eos bx - eos ax 


5 1 


a 2 -¿> 2 


2 +a 2 )(s 2 +6 2 ) 


51. 


a sin ax — Z> sin bx 
a 2 -b 2 


„2 
5 


(s 2 +a 2 )(s 2 +b 2 ) 


52. 


a cosax -b~ cosbx 

a 2 -b 2 


* 3 


(s 2 +a 2 )(s 2 +b 2 ) 


53. 


b sinh ax - a sinh bx 
a 2 -b 2 


ab 


(s 2 -a 2 ){ S 2 -b 2 ) 


54. 


cosh ax - cosh bx 
a 2 -b 2 


s 


(s 2 -a 2 )(s 2 -b 2 ) 
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A*) 


F(s) = £{Kx)) 


55. 


a sinh ax-b sinh bx 
a 2 -b 2 


s 2 


(s 2 -a 2 )(s 2 -b 2 ) 


56. 


a 2 cosh ax-b 2 cosh/).* 
a 2 -b 2 


s 2 


(s 2 -a 2 )(s 2 -b 2 ) 


57. 


1 . 
x sin ax 

a 


a 2 


s 2 (s 2 +a 2 ) 


58. 


— sinh ax - x 
a 


a 2 


s 2 (s 2 -a 2 ) 


59. 


, ax . 
1 - eos ax sin ax 

2 


a 4 


s(s 2 +a 2 ) 2 


60. 


1 - cosh ax + — sinh ax 
2 


a 4 


s{s 2 -a 2 ) 2 


61. 


. b 2 cosax-a 2 cosbx 

I+ „>-*> 


a 2 b 2 


s(s 2 +a 2 )(s 2 +b 2 ) 


62. 


1 + 


b cosh ax - a cosh bx 


a 2 b 2 




a ¿ -b ¿ 


s(s 2 -a 2 )(s 2 -b 2 ) 


63. 


1 
8 


(}>-a 2 x 2 )smax - 
3ax eos ax 




a 5 


{s 2 +a 2 f 


64. 


— sin ax - ax eos ax 

8 L J 


a s 


(s 2 +a 2 ) 3 


65. 


1 

8 


r ?7* 

(l + ú! A" )sin&x - 
ax eos ax 




a 3 s 2 


(s 2 +a 2 ) 3 
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M 


ffe) = *{/(*)} 


66. 


i 

8 


(3 + a 2 x 2 ) sinh ax - 
3ax cosh ax 




a 5 


(* 2 -* 2 ) 3 


67. 


X 

— (ax cosh ax - sinh ax) 
8 


a 3 j 


(* 2 -* 2 ) 3 


68. 


1 

8 


axcoshax - 

(1 -<2 2 x")sinhíix 




a 3 , 2 


is 2 -a 2 ? 


69. 


i-a-e -1 '')" 


1 


s (a5 + \)(as + 2)---(as + n) 


70. 


sin (ax + b) 


ssiab + a eos 6 

s 2 +a 2 


71. 


eos (a* + b) 


scosb-asinb 

2 2 

s +a 


72. 


e -ax_ e axf2 


ax*j3 
eos 

2 

/r . ax-JÍ 
V3 sin 

V 2 , 


/ 


3a 2 
s 3 + a 3 


73. 


\ + 2ax 

■\J7DC 


s + a 


sjs 


74. 


é^/Jiñ 


1 


V-s + a 


75. 


1 (e bx e") 

IXyJTDC 




yjs-a -*Js—b 


76. 


1 ^r~ 
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A*) 


F(s) = £{f{x)} 


77. 


, — cosh 2-Jax 

->J7ZX 


— e als 

4i 


78. 


1 • 1¡— 

—j^=sm2«Jax 


-3/2 -a/s 

s e 


79. 


. — sinh 2-y/ax 


-3/2 a/s 

s e 


80. 


J Q {2-Jax) 


-e~ a/s 
s 


81. 


yjxf aJ l (2jax) 


— e~ a/s 
s 2 


82. 


(x!a) {p ' x)l2 J p _ l {2^x') 
(p>0) 


—p —a/s 

s e 


83. 


Jo (x) 


1 


V* 2 +l 


84. 


A(x) 


js 2 +l-s 


V* 2 +l 


85. 


J p (x) (p > -1) 


(Vj 2 +1-í) p 


V* 2 +i 


86. 


x p J p {ax) \p>-~A 


(2a)T( P+ i) 


V^ 2 +*y +a/2) 


87. 


x p ~ l 

=- (P>0) 

r(p) 


i 

s p 


88. 


4 "" ! Y i»-0/2) 

(2n)¡V^ 


1 


5 W V7 
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A*) 


F(s) = £{f{x)j 


89. 


x p ~ ] 




1 




rf e (P>0) 

r(p) 


(s + a) p 




90. 


1-e" 

X 


s-a 

ln 

s 


91. 


bx „ax 

e -e 


ln 5 "* 
s-b 


X 


92. 


— sinh ax 

X 


, 5 + a 
ln 


93. 


2 

— (l-cosax) 

X 


. s 2 +a 2 


94. 


— (cosbx-cosax) 

X 


, s 2 +a 2 


95. 


sin a* 


a 

arctan — 

s 


X 


96. 


2 ■ u 
— sin ax eos bx 

X 


las 
arctan — — 

s 2 -a 2 +b 2 


97. 


sin \ax\ 




a 


í 1 + e -(^a )S \ 
\ L e ) 




2 2 
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Equations (cont.) 

n-th order, 34-35 

ordinary differential, 2 

partial differential, 2-3 

second-order, 34, 35-36 

separable, 7 
Equilibrium position, 47 
Euler's method, 106, 108, 111 
Exact equations, 7, 10-12 

Falling body problems, 21-22 
First-order differential equations, 
8-14 

applications, 20-26 

numerical methods, 105-112 
First-order systems, 79-83, 110 
Forced motion, 53 
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first-order, 105-112 
homogeneous, 30, 34-37, 39- 

44 
nonhomogeneous, 39-44 
reduction to first-order system, 

79-83 
second-order, 47-53, 86, 115- 

117 
solutions, 34-37, 79-83 
theory of solutions, 29-32 
Linearly dependent solutions, 3 1 
Linearly independent solutions, 

31 
Linear systems, 67 
Line elements, 1 05 

Matrices, 70-75 
Matrix addition, 71 
Matrix exponential, 73-75 
Matrix multiplication, 71 
Matrix solution methods, 79-83 
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Method of Frobenius, 90-91 
Method of undetermined coeffi- 

cients, 40-42 
Methods 
Adams-Bashforth-Moulton, 

109,111-112 
Euler's, 106, 108, 111 
Frobenius, 90-91 
modified Euler's, 108 
numerical for first-order equa- 

tions, 105-112 
predi ctor-corrector, 108 
Runge-Kutta, 108-109, 111 
undetermined coefficients, 40- 
42 
Modified Euler's method, 108 
Motion, 52-53 

Natural frequency, 53 
Natural length, 49 
Newton's law of cooling, 21 
Newton's second law of motion, 

22 
Nonhomogeneous equations, 32, 

39-44, 88-89 
Notation, 3 

nth-order equations, 34-36 
Numerators, 59 
Numerical instability, 107 
Numerical methods, 105-112 

Ordinary differential equations, 2 
Ordinary points, 87 
Orthogonal trajectories, 26 
Oscillatory damped motion, 53 
Overdamped motion, 53 

Partial differential equation, 2-3 
Partial fractions, 58-59 



Particular solutions, 4 

Period, 53 

Phase angle, 53 

Power series, 86-92 

Power series method, 87- 88 

Predi ctor-corrector method, 108 

Problems 

boundary- valué, 5, 115-120 

buoyancy, 51-52 

decay, 20-21 

dilution, 23-25 

eigenvalue, 117 

electrical circuit, 49-51 

falling body, 21-22 

growth, 20-21 

initial- valué, 44, 82 

spring, 47-49 

Sturm-Liouville, 117-118 

temperature, 21 

Recurrence formula, 88 
Regular singular points, 89-90 
Runge-Kutta method, 108-109, 
111 

Scalar multiplication, 71 
Second-order equations, 34, 35- 

36 
Second-order linear differential 
equations, 47-53, 86, 115- 
117 
Separable equations, 7 
Simple harmonic motion, 53 
Solutions, 3-5 
classifying, 52-53 
first-order differential equa- 
tions, 8-14 
initial- valué problems, 82 
Laplace transforms, 65-67 



136 DIFFERENTIAL EQUATIONS 



Solutions (cont.) 
linear differential, 34-37, 79- 

83 
linearly dependent, 3 1 
linearly independent, 3 1 
linear systems, 67 
matrix methods, 79-83 
no initial conditions, 82-83 
particular, 4 
power series, 86-92 
theory of, 29-32 

Spring problems, 47-49 

Square matrix, 70 

Standard forms, 5-6 

Steady-state motion, 53 

Step sizes, 107 

Sturm-Liouville problems, 117- 
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Temperature problems, 21 
Transient motion, 53 
Translations, 61-62 

Undamped motion, 53 
Underdamped motion, 53 
Undetermined coefficients, 40-42 
Unit step function, 60-61 

Variable coefficients, 30, 86 
Variation of parameters, 42-44 
Vectors, 70 

Wronskian, the, 31-32 



